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Abstract Largely motivated by the development of highly sensitive gravitational-wave detec- 
tors, our understanding of merging compact binaries and the gravitational waves they generate 
has improved dramatically in recent years. Breakthroughs in numerical relativity now allow us 
to model the coalescence of two black holes with no approximations or simplifications. There 
has also been outstanding progress in our analytical understanding of binaries. We review these 
developments, examining merging binaries using black hole perturbation theory, post-Newtonian 
expansions, and direct numerical integration of the field equations. We summarize these ap- 
proaches and what they have taught us about gravitational waves from compact binaries. We 
place these results in the context of gravitational-wave generating systems, analyzing the impact 
gravitational wave emission has on their sources, as well as what we can learn about them from 
direct gravitational-wave measurements. 
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1 Introduction 

1.1 History and motivation 

Most physics students learn to solve the Newtonian gravitational two-body prob- 
lem early in their studies. An exact solution for point masses requires only a few 
equations and an elliptic integral. Coupling this simple solution to perturbation 
theory lets us include the effect of non-spherical mass distributions and the im- 
pact of additional bodies. We can then accurately model an enormous range of 
astrophysically interesting and important systems. 

By contrast, no exact analytic solution describes binaries in general relativ- 
ity (GR). GR is nonlinear (making analytic solutions difficult to find except for 
problems that are highly symmetric) and includes radiation (so that any bound 
solution will evolve as waves carry off energy and angular momentum). Indeed, 
for systems containing black holes, GR doesn't so much have a "two-body" prob- 
lem as it has a "one-spacetime" problem: One cannot even delineate where the 
boundaries (the event horizons) of the "bodies" are until the entire spacetime 
describing the binary has been constructed. 

Such difficulties in describing binary systems in GR bedeviled the theoretical 
development of this topic. Many early discussions centered on the even more 
fundamental question of which motions would generate radiation and which would 
not. A particularly clear formulation of the confusion is expressed in attempts to 
answer the following question: // a charge falls freely in the Earth's gravitational 
field, does it radiate? On one hand, in this non-relativistic limit, we should expect 

3 
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that our usual intuition regarding accelerati ng charges wou ld hold, and a falling 



JacksonI (jl975l ) with an acceleration 



charge should radiate exactly as described in 
a = —gcr- On the other hand, in GR a falling charge simply follows a geodesic 
of the Earth's curved spacetime; it is not "accelerating" relative to freely falling 
frames, and so is not accelerating in a meaningful sense. The charge just follows 
the trajectory geometry demands it follows. In this picture, the falling charge 
appears to not radiate. John Wheeler once asked a group of relativity theorists 
to vote on whether the falling cha rge radiates or n ot; their responses were split 



almost precisely down the middle (jKennefick 



20071 . p. 157jl| 



Similar conceptual issues affect the 



general two-body p roblem in relativity. As 



recently as 1976, it was pointed out (jEhlers et al 



19761 ) that there had not yet 



been a fully self-consistent derivation for the energy loss from a binary due to 
gravitational-wave (GW) backreaction. A major mo tivation for this criti cism 



1975|). It 



was the discovery of the binary pulsar PSR 1913+16 (|Hulse Tavloij 
was clear that this system would be a powerful laboratory for testing theories 
of gravity, including the role of GW emission. However, as Ehlers et al. spelled 
out, the theoretical framework needed for such tests was not in good shape. 
Various approaches to understanding the evolution of binary systems tended to 
be inconsistent in the nature of the approximations they made, often treating the 
motion of the members of the binary at a different level of rigor than they treated 
the solution for the spacetime. These discrepancies weie most notable when the 



Dewitt fc Brehmg l| 196(3 ) : The charge does 



^The correct answer is now understood thanks to 
radiate, precisely reproducing the non-relativistic hmit. The intuition that the charge follows a 
geodesic is not quite right. Though the charge "wants" to follow a geodesic, the charge plus its 
associated field is extended and nonlocal, and so cannot follow a geodesic. The bending of the 
charge's field as it falls in spacetime enforces the laws of radiation emission. 
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members of the binary were strongly self gravitating; a weak-field approach is 

ill-suited to a binary containing neutron stars or black holes. These calculations 

generally predicted that the system would, at leading order, lose energy at a 

rate related to the third time derivative of the source's "quadrupole moment." 

However, the precise definition of this moment for strong- field sources was not 

always clear. 

The Ehlers et al. criticism served as a wake-up call, motivating the formu- 
lation of methods for modeling binary dynamics in a mathematically rigorous 
and self-consistent manner. Several independent approaches were developed; a 
concise and coge nt summary of the principles and the leading results is given by 



Damourl (|l983bl ). For the purpose of this review, a major lesson of the theoretical 
developments from this body of work is that the so-called "quadrupole formula" 
describing the loss of orbital energy and angular momentum due to GW emission 
is valid. Somewhat amazingly, one finds that the equations of motion are insen- 
sitive to much of the detailed structure of a binary's members. Features such as 
the members' size and compactness ca n be absorbed i nto the definition of their 



masses. This "principle of effacement" (jP amour 



1983al ) tells us that the motions 



of two bodies with masses rrii and m2 can be predicted independent of whether 
those bodies are stars, neutron stars, or black holeqj. Over 30 years of study have 
since found extraordinary agreem ent between prediction an d observation for the 



evolution of PSR 1913-M6's orbit (jWeisberg k Tavlor 



20051 ) . Additional inspiral- 



ing systems have been discovered; in all cases for which we have enough data to 



^Other aspects of the members' structure cannot be so simply absorbed by the principle of 
effacement. For example, at a certain order, the spins of a binary's members impact its motion. 
Spin's effect s cannot be absorbed int o the definition of mass, but affect the binary's dynamics 
directly. See 



Thorne fc Ifartld l|l985r ) for further detailed discussion. 
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discern period evolution, the data agree wi t h theory to within r neasurement pre 



cision 



2008 



Stairs et al 



1998 



Nice et al 



Bhat. Bailes &: Verbiest 



2005 



Jacobv et al 



2006 . 



Kramer &: Stairs 



20081 ). At least one additional rece ntly discovere d 



2008|) 



system is likely to show a measurable inspiral in the next few years (jKasianl l! 
These measurements have validated our theory of GW generation, and are among 
our most powerful tests of GR. 

Measuring GWs with the new generations of detectors will require even more 
complete models for their waveforms, and hence complete models of a binary's 
evolution. Motivated by this, our theoretical understanding of merging compact 
binaries and their GWs has grown tremendously. The purpose of this review 
is to summarize what we have learned about these systems and their waves, 
focusing on theory but connecting it to current and planned observations. We 
examine the analytic and numerical toolkits that have been developed to model 
these binaries, discuss the different regimes in which these tools can be used, and 
summarize what has been learned about their evolution and waves. We begin by 
first examining compact binaries as astrophysical objects, turning next to how 
they are treated within the theory of GR. 



1.2 Compact binaries: The astrophysical view 

From the standpoint of GR, all compact binary systems are largely the same 
until their members come into contact, modulo the value of parameters such as 
the members' masses and spins and the orbital period at a given moment. This 
means that we only need one theoretical framework to model any compact binary 
that we encounter in nature. From the standpoint of astrophysics, though, all 
compact binary systems are not the same: A 1.4 Mq — 1.4 Mq neutron star binary 
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forms in very different processes than those which form a 10^ Mq — 10^ Mq black 

hole binary. In this section, we summarize the astrophysical wisdom regarding 

the various compact binary systems that should be strong generators of GWs. 

Compact binaries are organized most naturally by their masses. At the low end 

we have stellar-mass binaries, which include the binary pulsars discussed in the 

previous section. The data on binaries in this end are quite solid, thanks to the 

ability to tie models for the birth and evolution of these systems to observations. 

At least some fraction of short gamma-ray bursts are likely to be associated with 



the mergers of ne utron star-neutron star (NS-IN 



(BH-NS) systems (|Eichler et al.l 



19891 : 



Fox et al 



S) or black hole-neutron star 



20051 ). Gamma-ray telescopes 



may already be telling us ab out compact binary merger many times per year 



^Nakar. Gal- Yam k Fox 



200a). 

There is also evidence that nature produces supermassive binaries, in which 
the members are black holes with M ~ 10^ — 10^ Mq such as are found at 
the centers of galaxies. As described in more detail below, theoretical argu- 
ments combining hierarchical galaxy growth scenarios with the hypothesis that 
most galaxies host black holes generically predict the formation of such bina- 
ries. We have now identified quite a few systems with properties indicating that 
they may host such binaries. The evidence inc l udes active galaxies wi th double 



cores ( Komossa et al. 



20031 . 



Maness et al 



with doubly-peaked emission lines ( 



2004 



Zhou et al. 



Rodriguez et al 



2004 



2006ll 



Gerke et al 



systems 



2003); hehcal 



Beeelman. Blandford & Rees 


1980. 


Conwav & Wrobel 


1995 


Lobanov & Roland 



20051 ) : and systems that 



OJ287 (Valtonen et al 



appe ar to be periodic or semi-periodic, such as the blazar 



20081 ). There are also sources which suggest the system 
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hosted a binary that recently merged, leading to the spin flip of a radio jet 



(Merritt fc Ekers 



2002) or to the interruption and later restarting of accretion 



activity (|Liu. Wu &: Cao 



20031 ) ■ As surveys go deeper and resolution improves, 
we may expect the catalog of candidate supermassive black hole binaries to ex- 
pand. 

Turn now from the observational evidence to theoretical models. If we as- 
sume that our galaxy is typical, and that the inferred density of NS-NS systems 
in the Milky Way should carry over to similar galaxies (correcting for factors 
such as typical stellar age and the proportion of stars that form neutron stars), 
then we can estimate th e rate at w hich binary syst ems merge in the universe. 



Naravan. Piran &: Shemil (|l99ll ) and 



PhinnevI (|l99ll ) first performed such esti- 



mates, finding a "middle-of-the-road" estimate that about 3 binaries per year 
merge to a distance of 200 Mpc. More recent calculations based on later surveys 
and observations of NS-NS systems have amended this number somewhat; the 
total number expected to be measured by advanced detectors is in the range of 



several tens per y ear (see, e.g., 



Kalogera et al 



(120071') for a detailed discussion of 



1 

in presd ) for a summary). 



methodology, and iKim. Kalogera &: Lorimeii (ji 

Population synthesis gives us a second way to model the distribution of stellar 
mass compact binaries. Such calculations combine data on the observed distri- 
bution of stellar binaries with models for how stars evolve. Begin with a pair of 
main sequence stars. The more massive star evolves into a giant, transfers mass 
onto its companion, and then goes supernova, leaving a neutron star or black 
hole. After some time, the companion also evolves into a giant and transfers 
mass onto its compact companion (and may be observable as a high-mass x-ray 
binary). In almost all cases, the compact companion is swallowed by the enve- 
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lope of the giant star, continuing to orbit the giant's core. The orbiting compact 

body can then unbind the envelope, leaving the giant's core behind to undergo a 



super nova explosion and form a compact remnant. See 



Tauris Sz van den Heuvel 



(|2006l ). especially Fig. 16.12, for further discussion. 

An advantage of population synthesis is that one can estimate the rate of 
formation and merger for systems which we cannot at present observe, such as 
stellar mass black hole-black hole (BH-BH) binaries, or for which we have only 
circumstantial evidence, such as neutron star-black hole (NS-BH) binaries (which 
presumably form some fraction of short gamma ray bursts). A disadvantage is 
that the models of stellar evolution in binaries have many uncertainties. There 
are multiple branch points in the scenario described, such as whether the binary 
remains bound following each supernova, and whether the binary survives com- 
mon envelope evolution. As a consequence, the predictions of calculations based 
on population synthesis can be quite diverse. Though different groups generally 
agree well with the rates for NS-NS systems (by design), their predictions for 



NS-B 



1997 



i and BH-BH systems diffe r by quite a bit (jYungelson &: Portegies Zwart 



Postnov &: Prokhorov 



19991 ). New data are needed to clear the theoretical 
cobwebs. 

Binaries can also form dynamically through many-body interactions in dense 
environments, such as globular clusters. In such a. cluster, the most massive 



19691 ): as such, the 



bodies will tend to sink through mass segregation (jSpitzeii 
core of the cluster will become populated with the heaviest bodies, either stars 
which will evolve into compact objects, or the compact objects themselves. As 
those objects interact with one another, they will tend to form massive bi- 
naries; calculations show that the production of BH-BH binaries is particu- 
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larly favored. It is thus likely that globula r clusters will act as "engines" for 



the production of massive compact binaries (IPortegies Zwart &: McMillan 



Q'Learv. O'Shaughnessv &: Rasio 



20071 . 



Mackev et al 



2000, 



200S ). 



As mentioned above, the hierarchical growth scenario for gal axies, coupled with 



the hypothesis that m ost galactic bulges host large black holes ([Kormendv &: Gebhardt 



2001 



Ferrarese 



2003) generically predicts the formation of supermassive binaries, 



especially at high redshifts when mergers were common. The fi rst careful discus 



sion of this scenario was bv lBegelman. Blandford &: Reed (|l98d ). In order for the 



black holes to get close enough to merge with one another due to GW emission, 
the black holes hosted by the merging galaxies must sink, via dynamical friction, 
into the center of the newly merged object. The binary thus formed will typically 
be very widely separated, and only harden through interactions with field stars 
(ejecting them from the center). For some time, it was unclear whether there 
would be enough stars to bring the holes close enough that they would be strong 
GW emitters. It is now thought that, at least on the low end of the black hole 
mass function (M ^ 10^ — IO'^Mq), this so-called "last parsec problem" is not 
such a problem. Quite a few mechanisms have be en found to carry the binary 



into the regime where GWs c an merge the binary (lArmitage &: Natarajan 



Merritt &: Milosavljevic 



2002 



20051 ) . It is now fairly common to assume that some 



mechanism will bring a binary's members into the regime where they can merge. 



blac 


<. holes and galaxies in 


hierarchical scenarios ( 


Menou, Haiman &; Naravanan 


2001 


, 


Yu & Tremaine 


2002 




Volonteri. Haardt & Madau 


2003 


). Galaxy mergers 



appear to be a natural mechanism to bring "fuel" to one or both black holes, 
igniting quasar activity; the formation of a binary may thus be associated with 
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the duty cycle of qua sars (jHaiman. Ciotti &: Ostrikei 



Di Matteo et al 



2004 



Hopkins et al 



2005, 



20081 ) . Such scenarios typically find that most black hole mergers 
come at fairly high redshift (z ^ 3 or so), and that the bulk of a given black hole's 
mass is due to gas it has accreted over its growth. 

A subset of binaries in the supermassive range are of particular interest to 
the relativity theorist. These binaries form by the capture of a "small" (1 — 
100 M0) compact object onto an orbit around the black hole in a galactic cen- 
ter. Such binaries forni dynamically through stellar interactions in the core 



(jSigurdsson &: Rees 



19971 . 



Hopman &: Alexander 



2OO5I ): the formation rate pre- 



dicted by most models is typical 



per year (jHopman &: Alexander 



10 ^ extreme mass ratio binaries per galaxy 



2OO5I ). If the inspiraling object is a white dwarf 



producing an x-ray or gamma- 


ray flare ( 


Rathore. Blandford &; Broderick 


2005, 


Menou. Haiman & Kocsis 


2008 




Rosswoe. Ramirez-Ruiz Sz Hix 


2008 


). If the in- 



spiraling object is a neutron star or black hole, it will be swallowed whole by the 
large black hole. As such, it will almost certainly be electromagnetically quiet; 
however, as will be discussed in more detail below and in Sec. HI its GW signature 
will be loud, and is a particularly interesting target for GW observers. 



1.3 Compact binaries: The relativity view 

Despite the diverse astrophysical paths to forming a compact binary, the end 
result always looks more-or-less the same from the standpoint of gravity. We 
now briefly outline the general features of binary evolution in GR. As described 
near the beginning of Sec. [H in GR a binary is not so much described by "two 
bodies" as by "one spacetime." The methods used to describe this spacetime 
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depend on the extent to which the two-body description is useful. 

Although it is something of an oversimplification, it is use ful to divide the evo- 
l ution of a binary into two or three broad epochs, following 



Flanagan &: Hughes 



(|l998l l. First is the binary's inspiral, in which its members are widely separated 
and can be readily defined as a pair of distinct bodies. During the inspiral, 
the binary's mean separation gradually decreases due to the secular evolution of 
its orbital energy and angular momentum by the backreaction of gravitational 
radiation. The bodies eventually come together, merging into a single highly 
dynamical and asymmetric object. We call the merger the final transition from 
two bodies into one. If the final state of the system is a black hole, then its last 
dynamics are given by a ringdown as that black hole settles down from the highly 



distorted post-m erger state into a Kerr 



theorems of GR ( Carter 



1971 



Robinson 



jlack hole as required by the "no hair" 



1973 ). 



How we solve the equations of GR to model a binary and its GWs depends 
upon the epoch that describes it. When the binary is widely separated, the post- 
Newtonian (pN) expansion of GR works very well. In this case, the Newtonian 
potential cp = GM/rc^ (where M = mi + m2 is the total mass of a binary, and r 
is the orbital separation) can be taken to be a small parameter. Indeed, we must 
always have r ^ (a few) x GM/c^: The closest the members of the binary can 
come is set by their physical size, which has a lower bound given by the radius 
they would have if they were black holes. The pN expansion is what we get by 
iterating GR's field equations from the Newtonian limit to higher order in cf). We 
review the basic principles of the pN expansion and summarize important results 
in Sec. \5\ 

Some binaries, such as the extreme mass ratio captures described in Sec. 11.21 
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will have nii <^ m2- For these cases, the reduced mass ratio rj = fi/M = 

mim2/{mi + 7712)'^ can define a perturbative expansion. In this limit, one can 
treat the binary's spacetime as an exact background solution with mass M (e.g., 
that of a black hole) perturbed by a smaller mass ^u. By expanding GR's field 
equations around the exact background, one can typically derive tractable equa- 
tions describing the perturbation and its evolution; those perturbations encode 
the dynamical evolution of the binary and its evolution. We discuss perturbative 
approaches to binary modeling in Sec. [H 

For some binaries, no approximation scheme is useful. Consider, for example, 
the last moments of two black holes coming together and fusing into a single body. 
In these moments, the spacetime can be highly dynamical and asymmetric; no 
obvious small parameter organizes our thinking about the spacetime. Instead, one 
must simply solve Einstein's field equations as exactly as possible using numerical 
methods. The essential question this field of numerical relativity asks is how 
one can take a "slice" of spacetime (that is, a single 3-dimensional moment of 
time) and use the field equations to understand how that spacetime evolves into 
the future. This requires us to explicitly split "spacetime" into "space" and 
"time." Progress in numerical relativity has exploded in recent years. Since 
2005, practitioners have moved from being able to just barely model a single 
binary orbit for a highly symmetric system to modeling multiple orbits and the 
final coalescence for nearly arbitrary binary masses and spins. We summarize the 
major principles of this field and review the explosion of recent activity in Sec.[6j 

Roughly speaking, for comparable mass binaries, pN methods describe the in- 
spiral, and numerical relativity describes the merger. The line dividing these 
regimes is fuzzy. A technique called the effective one-body (FOB) approximation 
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blurs it even further by making it pos sible to extend pN techniques beyond their 



naive range of validity (jDamour 



|2008| ). at least when both members of the binary 
are black holes. [When at least one of the members is a neutron star, at some 
point the nature of the neutron star fluid will have an impa c t. Det a iled numerica. 



model ii ig will then surely be critical; se e 



mm 



Baiotti. Giacomazzo fc Rezzollal (|2008l ). and 



Shibata &: UrvijI (1200611 . 



Shibata et al 



Itienne et al 



2009|) for ex- 



amples of recent progress.] Detailed tests show that using EOB methods greatly 
extends the domain for which ai ialytical waveform mo d els can accurately ru odel 



numerical relativity waveforms ([Buonanno et al 



2007 



D amour et al 



2003). A 



brief discussion of these techniques is included in Sees. [5] and [H 

Finally, it's worth noting that the ringdown waves that come from the last 
dynamics of a newly-born black hole can also be modeled using perturbation the- 
ory. The spacetime is accurately modeled as a black hole plus a small deviation. 
Perturbation theory teaches us that black holes "ring" in a series of modes, with 
frequencies and damping tim es that are controlled by the mass and the spin of 



the black hole (Leaver 



19851 ). Any deviation from an exact black hole so^ 



is carr ied away by such modes, enforcing the black hole no-hair theorems (jPrice 



1972a 



ution 



b|). We will not say much about the ringdown in this review, except to 
note that the last waves which come from numerical relativity simulations have 
been shown to agree excellently with these perturbative calculations. 



1.4 Notation and conventions 

The underlying theory of GWs is general relativity (GR); we review its most 
crucial concepts in Sec. [2j Because multiple conventions are often used in the 
GR literature, we first describe our notation and conventions. 
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Throughout this review, Greek indices denote spacetime components of tensors 

and vectors. These indices run over (0, 1, 2, 3), with generally denoting time t, 

and (1, 2, 3) denoting spatial directions. Spacetime vectors are sometimes written 

with an overarrow: 

A = {A^'} = {A°, A\A^,A^) . (1) 

Equation ([1]) should be read as "The vector A has components whose values in 
a specified coordinate system are A^, A^, A^, and A'^." Lowercase Latin indices 
denote spatial components. Spatial vectors are written boldface: 

SL = {a'} = {a\a'^,a^) . (2) 

We use the Einstein summation convention throughout, meaning that repeated 
adjacent indices in superscript and subscript positions are to be summed: 

3 

A'^B^^Y.A^B^. (3) 

Indices are raised and lowered using the metric of spacetime (discussed in more 
detail in Sec. [2]) as a raising or lowering operator: 

A^^B^ = g^.A^B'' = A,B^ = g^'A^B, . (4) 

When we discuss the linearized limit of GR (particularly in Sec. [3]), it is useful 
to work in coordinates such that the spacetime metric can be written as that of 
special relativity plus a small perturbation: 

where rj^^, = diag( — 1, 1, 1, 1). This means that the spatial part of the background 
metric is the Kronecker delta 6ij. For certain calculations in linearized theory, 
it is useful to abuse the Einstein summation convention and sum over repeated 



16 Scott A. Hughes 

adjacent spatial indices regardless of position: 

3 

aj6* = aj6* = Oibi = a^b^ . (6) 

i=l 

This is allowable because using the Kronecker delta for the spatial part of the 
metric means a* = to this order. 

Throughout this review, we abbreviate the partial derivative 

d 



■ (7) 



With this notation defined, we can write = g^'^du. 

Finally, it is common in GR research to use units in which the gravitational 
constant G and the speed of light c are set to 1. This has the advantage that 
mass, space, and time have the same units, but can be confusing when applied to 
astrophysical systems. In deference to the astronomical audience of this review, 



we have put Gs and cs back i nto the relativity formu ^ 



that G = c = 1 everywhere is 



as. A n exception to the rule 



Shapiro fc Teukolskvl (jl983l ). especially Chap. 15. 



2 Synopsis of general relativity 

GR describes gravity as geometry. The foundation of this is the metric, which 
provides a notion of spacetime distance. Suppose event A occurs at coordinate 
x", and event B at + dx°'. The proper spacetime separation, ds, between A 
and B is given by 

ds"^ = ga^dx'^dx^ . (8) 

The metric gaf3 translates the information in coordinates, which can be arbitrary, 
into a "proper" quantity, which can be measured. In the limit of special relativity, 
ga/3 becomes ijajj (defined in Sec. ll.4p . The general spacetime metric is determined 
by the distribution of mass and energy; we describe how to compute it below. It 
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will sometimes be useful to work with the inverse metric g'^^, defined by 

5"^ff/37 = 5% • (9) 

The metric also takes inner products between vectors and tensors: 

A-B = g^fsA''Bf^ . (10) 

A is timelike if ^ • ^4 < 0, spacelike if ^ • A > 0, and lightlike or null if A • ^4 = 0. 

Consider a worldline or spacetime trajectory z'^{t), where r is "proper time" 
(time as measured by an observer on that worldline). The vector u'^ = dz^/dr is 
the tangent to the worldline. If u'^ is timelike, it is the 4-velocity of an observer 
following the worldline, and is normalized u^U/j, = —1. Suppose the worldline 
extends from A to B. The total spacetime separation between these points is 

s = dT^go^pu^uP . (11) 

We now extremize s: fix the endpoints, allow quantities under the integral to 
vary, but require the variation to be stationary (so that 5s = Qi). The which 
extremizes s is given by the geodesic equation: 

^ + O^uV = 0. (12) 

We have introduced the "connection" F^^-y; it is built from the metric by 

r"/37 = ^^^'^ (^t5/./3 + dpg^^ - d^gp^) . (13) 

Geodesies are important for our discussion because freely falling bodies follow 
geodesies of spacetime in GR. Geodesies express the rule that "spacetime tells 
bodies how to move." 

Timelike geodesies describe material bodies. Null geodesies, for which u^^u^ = 
0, describe massless bodies or light rays. Our discussion above describes null 
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geodesies, with one modification: We cannot parameterize a null worldline with 

r, as proper time is not meaningful for a "speed of light" trajectory. Instead, 

one uses an affine parameter A which uniformly "ticks" along that trajectory. A 

convenient choice is to set = dx"/dX to be the 4-momentum of our radiation 

or massless body. With this choice, our discussion describes null trajectories just 

as well as timelike ones. 

The connection also defines the covariant derivative of a vector or tensor: 

VaA/S = dc^Ap-Ai^Vap , 
^a-^^j — doiA^j + A^^T^gifj^ — A^ nT^a-y . (14) 

The pattern continues as we add indices. This derivative follows by comparing 
vectors and tensors that are slightly separated by parallel transporting them to- 
gether to make the comparison; the connection encodes the twists of our curved 
geometry. Using the covariant derivative, the geodesic equation can be written 

n^V^jn" = . (15) 

In curved spacetime, nearby geodesies diverge. Because a geodesic describes a 
freely-falling body, the rate of at which geodesies diverge describes tides. Let 
be the displacement between geodesies. Then the rate of divergence is given by 

^ = R'^fs-ysnf'u^^' ■ (16) 
We have introduced the Riemann curvature tensor, 

R°'p^s = d^^'^ps — dsF'^fjj + T^'njT^ps — ^""nS^^p-y ■ (17) 
Some variants of Riemann are important. First, there is the Ricci curvature: 

Ra/3 = R^an/S ■ (18) 
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Ricci is the trace of Riemann. Taking a further trace gives us the Ricci scalar, 

= . (19) 
The Ricci tensor and Ricci scalar combine to produce the Einstein curvature: 

Gal3 = Ral3 — ^dapR ■ (20) 

The Riemann curvature satisfies the Bianchi identity, 

''^jRalSfii/ + ^ IsR-icxfiv + Vo-ii^g^^jy = . (21) 

By tracing over certain combinations of indices, the Bianchi identity implies 

V"G„;3 = , (22) 

a result sometimes called the "contracted" Bianchi identity. 

So far, we have mostly described the mathematics of curved geometry. We 
must also introduce tools to describe matter and fields. The most important tool 
is the stress-energy tensor: 

rpiiv _ -p\\xyi of momentum in the direction. (23) 

An observer who uses coordinates (i,^;*) to make local measurements interprets 
the components of this tensor as 

T** = Local energy density (24) 

T** = Local energy flux (times c) (25) 

T** = Local momentum density (times c) (26) 

T*-^ = Local momentum flux (times c^); T** acts as pressure. (27) 

[The factors of c in Eqs. (|25|) - (|27|) ensure that the components of T^'^ have the 
same dimension.] Local conservation of energy and momentum is expressed by 

V^T^"' = . (28) 
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In GR, we generally lose the notion of global energy conservation: We cannot 

integrate Eq. ()28p over an extended region to "add up" the total energy and mo- 
mentum. This is because V^T'^'^ is a spacetime vector, and in curved spacetime 
one cannot unambiguously compare widely separated vectors. 

Einstein's hypothesis is that stress energy is the source of spacetime curvature. 
If T^'^ is our source, then the curvature must likewise be divergence free. The 
contracted Bianchi identity (j22p shows us that the Einstein tensor is the curvature 
we need. This logic yields the Einstein field equation: 

G^u — Tfj^i/ . (29) 

The factor SirG/c^ guarantees that this equation reproduces Newtonian gravity 
in an appropriate limit. Note its value: 

^ = 8.26 X 10-50 ^^^^ . (30) 
c erg/cm 

It takes an enormous amount of energy density to produce spacetime curvature 

(measured in inverse length squared). Note that the reciprocal of this quantity, 

times c, has the dimensions of power: 

— = 3.63 X 10^^ ere/sec . (31) 

This is the scale for the power that is generated by a GW source. 

3 Gravitational-wave basics 

We now give a brief description of how gravitational waves arise in GR. Our 
purpose is to introduce the main ideas of this field, and also provide some results 
against which the more complete calculations we discuss later can be compared. 
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3.1 Leading waveform 

Begin with "weak" gravity, so that spacetime is nearly that of special relativity, 

= ValS + Kl3 ■ (32) 

Take the correction to flat spacetime to be small, so that we can linearize about 
r]ai3. Consider, for example, raising and lowering indices: 

h»P ^ g^^^g^^h^, = v^'^v'^'h^, + 0{h^) . (33) 

Because we only keep quantities to first order in h, we will consistently use the 
flat metric to raise and lower indices for quantities related to the geometry. 
Applying this logic repeatedly, we build the linearized Einstein tensor: 

Gcf^ = \ (c>„a^ V + Qpd^Ko^ - d^dph - UKp + r^^pUh - Vapd^'d'^h^,) , (34) 

where h = r]°'^hap is the trace of haf3, and □ = rj'^^dadp is the flat spacetime 
wave operator. 

Equation (j34p is rather messy. We clean it up in two steps. The first is pure 
sleight of hand: We introduce the trace-reversed metric perturbation 

ha/3 = ha/3 " ^Va/3h ■ (35) 

With this definition, Eq. (j39p becomes 

Ga^ = \ {dad^\p + dpd^h^a " nha/3 " r^a^d^&^h^^) . (36) 

Next, we take advantage of the gauge-freedom of linearized gravity. Recall that 
in electrodynamics, if one adjusts the potential by the gradient of a scalar, Ap, — > 
~ ^M^' f^^sn the field tensor F^y = dfj_Ay — dyA^j, is unchanged. In linearized 
GR, a similar operation follows by adjusting one's coordinates: If one changes 
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coordinates rc" — > + ^" (requiring 5^.^" 1), then 



Scott A. Hughes 



(37) 



One can easily show [see, e.g., 



Carroll! (|2004l ). Sec. 7.1] that changing gauge leaves 



the Riemann tensor (and thus all tensors derived from it) unchanged. 
We take advantage of our gauge freedom to choose so that 



3" V = 



(38) 



This is called "Lorenz gauge" in analogy with the electrodynamic Lorenz gauge 
condition d'^A^ = 0. This simplifies our Einstein tensor considerably, yielding 

1, 



Gal3 



-Dh 



(39) 



The Einstein equation for linearized gravity thus takes the simple form 



Uh 



levrG, 



■a/3 



(40) 



A ny linear equatio n of this form can be solved using a radiative Green's function 



[e.g. 



JacksonI (|l975l ). Sec. 12.11]. Doing so yields 



/lQ,/3(x,t) = — / ■ d X 



X — X' 



(41) 



In this equation, x is the "field point," where h^p is evaluated, and x' is a "source 
point," the coordinate that we integrate over the source's spatial extent. Notice 
that the solution at t depends on what happens to the source at retarded time 
t — |x — x'l/c. Information must causally propagate from x' to x. 

Equation (j4ip is formally an exact solution to the linearized Einstein field 
equation. However, it has a serious problem: It gives the impression that ev- 
ery component of the metric perturbation is radiative. This is an unfortunate 
consequence of our gauge. Just as one can choose a gauge such that an isolated 
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point charge has an oscillatory potential, the Lorenz gauge we have used makes 

all components of the metric appear radiative, even if they are statiqf . 

Fortunately, it is not difficult to see that only a subset of the metric represents 

the truly radiative degrees of freedom in all gauges. We will only quote the resu lt 



Flanagan Huehesl pOOSl l 



here; interested readers can find the full calculation in 
Sec. 2.2: Given a solution hap to the linearized Einstein field equations, only the 
spatial, transverse, and traceless components K^^ describe the spacetime's 
gravitational radiation in a gauge-invariant manner. (The other components 
can be regarded as "longitudinal" degrees of freedom, much like the Coulomb 
potential of electrodynamics.) Traceless means 

5i,hJ^ = ; (42) 

"transverse" means 

dih^^ = . (43) 

This condition tells us that hj^" is orthogonal to the direction of the wave's 
propagation. Expanding hj^ in Fourier modes shows that Eq. (|13]) requires hj^ 
to be orthogonal (in space) to each mode's wave vector k. 

Conditions (j42p and (j43p make it simple to construct Kf^ given some solution 
hij to the linearized field equations. Let Ui denote components of the unit vector 
along the wave's direction of propagation. The tensor 

Pij = Sij - UiUj (44) 



■^In the electromagnetic case, it is unambiguous which field components are radiative and 
which are static. Similarl y, one can a l ways t ell which Riemann curvature components are radia- 



tive and which are static. lEddingtorJ (|l922l ) appears to have been the first to use the curvature 



tensor to categorize the gravitational degrees of freedom in this way. 
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projects spatial components orthogonal to n. It is then simple to verify that 

hj^ = hki {PkiPi, - \PklP^^ (45) 

represents the transverse and traceless metric perturbation. 

The simplest example solution to Eq. (|45p can be built by going back to Eq. 
(j4ip and focusing on the spatial components: 

,_^.^^|M^^ld^.v. (40) 

Consider a very distant source, putting |x — x'| ~ i?: 

hj - ^ § / Tij (x', t - R/c)d?x' . (47) 

To proceed, we invoke an identity. Using the fact that V^T^u = goes to 
d^^T^y = in linearized theory and in our chosen coordinates, we have 

d'Ttt + d^Tjt = , d% + d% = . (48) 

Combine these identities with the fact that 3* = —dt; use integration by parts 
to convert volume integrals to surface integrals; discard those integrals by taking 
the surface outside our sources. We then find 

J T,,i^',t)d^x' = \^j x^'x^'Tu{^',t)d^x' ^ l^Ht) . (49) 

We have introduced here the quadrupole moment lij. This allows us to at last 
write the transverse-traceless waveform as 

It is straightforward to show that the trace / = la does not contribute to Eq. 
([50]) . so it is common to use the "reduced" quadrupole moment, 

Xij = lij — '^^ij^ ■ (51) 
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The waveform then takes the form in which it is usually presented, 

the quadrupole formula for GW emission. 

A more accurate approximation than |x — x'| ~ i? is 

|x - x'l ~ i? - n • x' , (53) 

where n is the unit vector pointing from x' to x. Revisiting the calculation, 
we find that Eq. (j52p is the first term in a multipola r expans i on fo r the radia- 



tion. Detailed formulae and notation can be found in IThornd (| 19801 ). Sec. IVA. 
Schematically, the resulting waveform can be written 

^ ^(A^d^ , B^d^l^^"" 

Rc^ dt^ c'+i dt^ J • ^ ' 

1=2 ^ ■' 

We have hidden various factorials of I in the coefficients Ai and Bi; the superscript 
"STF" means to symmetrize the result on any free indices and re move the trace . 
The symbol 2i stands for the Ith mass moment of the source. iThornd (|l98Cll ) 



precisely defines If, for our purposes, it is enough to note that it represents an 
integral of / powers of length over the source's mass density p: 

li^ j p{x'){x')U^x' . (55) 

The mass moment plays a role in gravitational radiation similar to that played 
by the electric charge moment in electrodynamics. The symbol Si describes the 
/th mass-current moment] it represents an integral of I powers of length over the 
source's mass-current J = pv, where v describes a source's internal motions: 

Si^ j p{x')v{x'){x'y d^x' . (56) 

Si plays a role similar to the magnetic moment. 
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The similarity of the multipolar expansion (j54p for GWs to that for electro- 
magnetic radiation should not be a surprise; after all, our linearized field equation 
()40p is very similar to Maxwell's equation for the potential Afj_ (modulo an extra 
index and a factor of four). One should be cautious about taking this analogy 
too far. Though electromagnetic intuition applied to linearized theory works well 
to compute the GWs from a source, one must go beyond this linearized form to 
understand deeper aspects of this theory. In particular, one must go to higher or- 
der in perturbation theory to see how energy and angular momentum are carried 
from a radiating source. We now sketch how this works in GR. 

3.2 Leading energy loss 

Electromagnetic radiation carries a flux of energy and momentum described by 
the Poynting vector, S = (c/47r)E x B. Likewise, electromagnetic fields generate 
stresses proportional to (|Ep + |Bp)/87r. The lesson to take from this is that 
the energy content of radiation should be quadratic in wave amplitude. Properly 
describing the energy content of radiation requires second-order perturbation 
theory. In this section, we will dis c uss th e key concepts and ideas in this analysis, 



which was first given by 



Isaacson! (|l968l ). 



Begin by writing the spacetime 

gal3 = 9ap + e/lQ:/3 + • (57) 

We have introduced a parameter e whose formal value is 1; we use it to gather 
terms that are of the same order. Note that now we do not restrict the background 
to be flat. This introduces a conceptual difficulty: Measurements directly probe 
only the total spacetime gap (or a derived surrogate like curvature), so how do we 
distinguish background gap from perturbation? The answer is to use separation of 
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lengthscales. Our background will only vary on "long" lengthscales and timescales 

C,T; our perturbation varies on "short" lengthscales and timescales A, r. We 

require £. ^ X and T ^ t. Let (/) denote a quantity / averaged over long scales; 

this averaging is well-defined even for tensors up to errors ©(A^/ZZ^). Then, to 

first order in e, 

9a/3 = {ga/s) (58) 
haP = gap - gal3 ■ (59) 

The second-order contribution will be of order /i^, and (as we'll see below) will 
have contributions on both long and short scales. 

Begin by expanding the Einstein tensor in e. The result can be written 

G^Aga/3) = G%{g^p) + eGl^iKp) + e\Gl,{Kp) + Gl,{j^p)] . (60) 

For simplicity, take the spacetime to be vacuum — there are no non-gravitational 
sources of stress and energy in the problem. We then require Einstein's equation, 
Gij,u = 0, to hold at each order. The zeroth order result, 

G%{gap) = , (61) 

is just a statement that we assume our background to be a vacuum solution. 
Expanding G^i/(/iq/3) = 0, we find 

- V - Rcf.pu'h^" = , (62) 

where □ = g'^'^'V^'Vu is the wave operator for g^i, (with the covariant deriva- 
tive on gfj.u), Rafj.f3u is the Riemann curvature built from cjf^^, and h^i, = h^y — 
{l/2)g^uh. Equation (j62p is just the wave equation for radiation propagating 
on a curved background. The coupling between h and the background Riemann 
tensor describes a correction to the "usual" geometric optics limit. 
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Next, consider second order: 

GUU) = -GlAKp) ■ (63) 

To make sense of this, invoke separation of scales. The second-order perturbation 
has contributions on both scales: 

= Ua/S) + Sjal3 , (64) 

where {ja/s) varies on long scales and Sja/s is oscillatory and varies on short scales. 
The second-order metric can now be written 

9al3 = 9al3 + + ^^^JaP , (65) 

where 

9al3=9af^ + (^'^{jaf^) , (66) 

is a "corrected" background which includes all pieces that vary on long scales. 

Now return to the Einstein equation (I60p . but consider its averaged value. 
Thanks to linearity, we can take the averaging inside the operator G^: 

{Gl,{Kp)) = Gl,{{Kp)) = Q. (67) 

We have used here {hap) = 0. Putting all this together, we find 

GlX9ap) + e\{Gl,{Kp)) + Gl,{{j^M = . (68) 

Let us rewrite this as 

G^.(9a/3 + e'Oa/3)) = -e^ (G^, (/i„^)) , (69) 

or, putting e = 1, 

GA9^0) = -{GUKp)) . (70) 
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Equation (j70p says that the second-order averaged Einstein tensor acts as a source 

for the long lengthscale background spacetime. This motivates the definition 

T^J^ = -^{GUf'ap)) ■ (71) 
Choosing a gauge so that Vf'h^^ = 0, r^w takes a very simple form: 



TS^ = ^{'^,Kp^.h"^)- (72) 



This quantity is known as the Isaacson stress-energy tensor ([Isaacson 



19681 ). 



3.3 The "Newtonian, quadrupole" waveform 

A useful exercise is to consider a binary with Newtonian orbital dynamics that 
radiates GWs according to Eq. (1521) . Further, allow the binary to slowly evolve 
by energy and angular momentum carried off in accordance with Eq. (j72p . 

Begin by considering such a binary with its members in circular orbit of sepa- 
ration R. This binary is characterized by orbital energy 

E = -miv^ + -m2V2 — = ^ , (73) 

(where M = mi + m2 and fi = mim2/M) and orbital frequency 



^orb = y-^. (74) 

Next consider the energy that GWs carry from the binary to distant observers. 
When evaluated far from a source, Eq. (|72|) gives a simple result for energy flux: 



dE _ 

dAdt ~ 327rG'^""*^ ^"'"'^ 



m^d^^^)r^ . (75) 



Plugging in Eq. ()52|) and integrating over the sphere, we find 

dt ~ J dAdt ~ 5c5 \ dt^ dt^ ■ ^ ' 
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For the Newtonian binary, 



1 
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(77) 



we choose coordinates for this binary such that the components of the separation 
vector are xi = RcosQorht, X2 = Rsmi^oj-\^t, X3 = 0. Inserting into Eq. (f76l) . we 
find 

(78) 



It ~ Tc^ 



We now assert that the binary evolves quasi-statically, meaning that any radiation 
carried off by GWs is accounted for by the evolution of its orbital energy: 



'dt ^~dt ~ ' 



(79) 



We evaluate dE°^^ /dt by allowing the orbital radius to slowly change in time, 

(80) 



dE_°'^ _ dE'^'^dR 
H ~ dR ~dt 



Combining Eqs. ([78]), ([79]), and ([80]), we find 



R{t) = 



256GVM^(tc-t) 
5c5 



1/4 



(81) 



This in turn tells us that the orbital frequency changes according to 



5c^ 



3/8 



(82) 



_256(GA^)5/3(t, -t)_ 

We have introduced the binary's chirp mass M = /i^/^M^/s, so called because it 
sets the rate at which the binary sweeps upward in frequency, or "chirps." We 
have also introduced the coalescence time tc, which formally describes when the 
separation goes to zero (equivalently, when the frequency goes to infinity). By 
rearranging Eq. (jSip . we find the time remaining for a circular binary of radius 
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R to coalesce due to GW emission: 

5 i?^ 

Tremaining = ^^^^ 



3 X lO'' years 
2 months 



2.8 Mq \ V -R ^ ^ 



M J \Rq 



M J \A\J 

„ / 2 X 10^ Mr..\^ f R , 
3 X lOVears — - . (83) 



M J VO-OOlpc 
The fiducial numbers we show here are for equal mass binaries, so // = M/4. 

Given the substantial eccentricity of many binaries, restriction to circular or- 
bits may not seem particularly realistic. Including eccentricity means that our 
binary will have two evolving parameters (semi- major axis a and eccentricity 
e) rather than just one (orbital radius R). To track their evoluti on, we must 



separately compute the radiated energy and angular momentum (jPeters 



Peters fc Mathews 



1963): 



1964 



-dt = B^^it^u^^ = y ^ ' 

dLz 2G .(PljjYi dPlkm\ 32 G n 4/->5 / n /or\ 

Because the binary is in the x — y plane, the angular momentum is purely along 
the z axis. The eccentricity corrections /(e) and g{e) are given by 

1 + 73 2 , 37 4 I , 7 2 

Using standard definitions relating the semi-major axis and eccentricity to the 



orbit's energy E and angular momentum L^, Eqs. (j84p and (j85p imply ([Peters 
1964 1 



dt 5 o' 

^ 304 



de 304 ^M^ 1 + ii^^ 



-e 



15 c5 a4 (1 - e2)5/2 ■ 



32 Scott A. Hughes 

It is then simple to compute the rate at which an eccentric binary's orbital period 

changes due to GW emission. This result is compared with data in investigations 
of GW generating binary pulsars. Because eccentricity tends to enhance a sys- 
tem's energy and angular momentum loss, the timescales given in Eq. (jSSp can 
be significant underestimates of a binary's true coalescence time. 

Notice that a binary's eccentricity decreases: GWs tend to circularize orbits. 
Many binaries are expected to be essentially circular by the time their waves 
enter the sensitive band of many GW detectors; the circular limit is thus quite 
useful. Exceptions are binaries which form, through capture processes, very close 
to merging. The extreme mass ratio inspirals discussed in Sec.|3]are a particularly 
interesting example of this. 

We conclude this section by writing the gravitational waveform predicted for 
quadrupole emission from the Newtonian, circular binary. Evaluating Eq. (j52|) . 
we find that hij has two polarizations. These are traditionally labeled "plus" and 
"cross" from the lines of force associated with their tidal stretch and squeeze. 
Taking our binary to be a distance D from Earth, its waveform is written 

h+ = ^ ^ (1 + COSn)cOs2$Ar(t) , 



= -^Tp- ( c3 J costsin2$Ar(t) , (89) 



where the phase 



^N{t) = j narhdt = $c 



c^t.-t) 



5/8 



(90) 



5GM 

and where / = {l/ir)d^]\f /dt is the GW frequency. The system's inclination l is 
just the projection of its orbital angular momentum, L, to the wave's direction 
of propagation n: cosi = L • n (where L = L/|L|). We show fiducial values for 
the GW amplitudes when we briefly describe GW measurement in Sec. [HI 
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In later sections, we will use Eq. (j89p as a reference to calibrate how effects we 

have neglected so far change the waves. Note that and depend on, and 

thus encode, the chirp mass, distance, the position on the sky (via the direction 

vector n), and the orientation of the binary's orbital plane (via L). 

3.4 Nonlinear description of waves 

The nonlinear nature of GR is one of its most important defining characteristics. 
By linearizing, perhaps we are throwing out the baby with the bathwater, failing 
to characterize important aspects of gravitational radiation. Fortunately, we can 
derive a wave equation that fu _ 
was apparently first derived by 



ly encodes all nonlinear features of GR. This 



Penrose ( 196 



RvanI (jl974l ) gives a very nice 



discussion of this equation's history and derivation. Begin by taking an additional 
derivative V''' of the Bianchi identity ()2ip , obtaining 

(91) 



where Ug = g^^V^Vs is a covariant wave operator. Next, use the fact that the 
commutator of covariant derivatives generates a Riemann: 



[V^, V^]^^ = (V^V^ - Vs^^)Pti = -R'^fi-fSPa 



[^jT^s]Pfiu — R pi-ySPcFV R i/^SP^cr • 



(92) 



(93) 



Extension to further indices is hopefully obvious. Manipulating (j9ip yields a 
wave equation for Riemann in w hich (Rieman n)^ acts as the source term. This 



is the Penrose wave equation; see 



RvanI (jl974l ) for more details. 



If spacetime is vacuum \T^i^ — 0, so that via the Einstein equation (|29p R^u 
0], the Penrose wave equation simplifies quite a bit, yielding 



^gRaf3fiu — '^RpLafirRv a "^R/iaar Ru /3 ~l~ R/iarcrR a(3 ■ (94) 



34 Scott A. Hughes 

A variant of Eq. ()94p underlies much of black hole perturbation theory, our next 

topic. 



4 Perturbation theory 

Perturbation theory is the first technique we will discuss for modeling strong-field 
merging compact binaries. The basic concept is to assume that the spacetime is an 
exact solution perturbed by a small orbiting body, and expand to first order in the 
binary's mass ratio. Some of the most interesting strong-field binaries have black 
holes, so we will focus on black hole perturbation theory. Perturbation theory 
analysis of binaries has two important applications. First, it can be a limiting 
case of the pN expansion: Perturbation theory for binaries with separations r ^ 
GM/c^ should give the same result as pN theory in the limit /x <^ M. We 
return to this point in Sec. O Second, perturbation theory is an ideal tool for 
extreme mass ratio captures, binaries created by the scattering of a stellar mass 
(m ~ 1 — 100 Mq) body onto a strong- field orbit of a massive {M ~ 10^ — lO'^ -^o) 
black hole. 

4.1 Basic concepts and overview of formalism 

At its most basic, black hole perturbation theory is developed much like the weak 
gravity limit described in Sec. 13. H replacing the flat spacetime metric r]a/3 with 
the spacetime of a black hole: 

9^.u = g^^ + Ku ■ (95) 

For astrophysical scenarios, one uses the Schwarzschild (non-rotating black hole) 
or Kerr (rotating) solutions for g^^. It is straightforward (though somewhat 
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tedious) to then develop the Einstein tensor for this spacetime, keeping terms 

only to first order in the perturbation h. 

This approach works very well when the background is non-rotating, 

(d,2)BH ^ gl^dx^'dx' = _(i_'^jAdt^^ , , + r^dn" , (96) 

V / (l-2M/rj 

where dO,"^ = dO^ + sin^ OdcfP' and M = GMjc? . We c onsider this special case in 



detail; our discussion is adapted from 



Rezzollal (j2003l ). Because the background 



is spherically symmetric, it is useful to decompose the perturbation into spherical 
harmonics. For example, under rotations in Q and 0, /iqo should transform as a 
scalar. We thus put 

/ioo = J^a/m(t,r)li„(0,0) . (97) 

Im 

The components /iQj transform like components of a 3- vector under rotations, and 
can be expanded in vector harmonics; hij can be expanded in tensor harmonics. 
One can decompose further with parity: Even harmonics acquire a factor ( — 1)' 
when (0, 0) [tt — 6,(p + tt); odd harmonics acquire a factor (— 1)'"*"-*^. 

By imposing these decompositions, choosing a particular gauge, and requiring 
that the spacetime satisfy the vacuum Einstein equation Gf^u = 0, we find an 
equation that governs the perturbations. Somewhat remarkably, the t and r 
dependence for all components of /i^jy for given spherical harmonic indices {l,m) 
can be constructed from a function Q{t,r) governed by the simple equation 

where r* = r + 2Af ln(r/2M — 1). The potential V{r) depends on whether we 
consider even or odd perturbations: 

2q{q + l)r^ + Gq^Mr'^ + ISqM^r + 18M^ 



^even(r) = (^1 - 



(qr + 3M 



(99) 
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where q = {I - + 2)/2; and 



, , / 2M 
Kdd(r) = 1 - — 



l{l + 1) 6M 



For even parity, Eq. ([98l) is the Zerilli equation (jZerillil ll97Cl ): for odd, it is 



the Reg ge- Wheeler equation (jRegge Wheeler 
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(100) 



19571 ). For further di scussion, 



i nclud ing how gauge is chosen and how to construct hnu from Q, see 



Rezzolla 



(|2003l ). Finally, note that when the spacetime perturbation is due to a body 
orbiting the black hole, these equations acquire a source term. One can construct 
the full solution for the waves from an orbiting body by using the source-free 
equation to build a Green's function, and then integrating over that source. 
How does this procedure fare for rotating holes? The background spacetime. 



{ds 



2\BU 



1 



2Mr 



dt^ 



AaMr sin^ 9 



P / P 

2 2 2Mra^sm^0^ 



6^ 

dtdcj) + ^dr^ + p^de^ 



+ + a' + 



P^ 



(101) 



where 



cM ' 



2 2,2 2 

p = r + a cos I 



A = - 2Mr + , 



(102) 



is now markedly nonspherical. [We have used "Boyer-Lindquist" coordinates, but 
the nonspherical nature is independent of coordinate choice. The spin parameter 
a has the dimension of mass, and must satisfy a < M in order for Eq. (jlOip 
to represent a black hole.] So, the decomposition into spherical harmonics is 
not useful. One could in principle simply expand G^,y = to first order in h^^ 
and obtain a partial differential equation in t, r, and 9. (The metric is axially 
symmetric, so we can easily separate the (j) dependence.) This author is unaware 
of such a formulatiory. One issue is that the gauge used for the iperturbatioii 



Since originally posting this paper, I was informed by Plamen Fiziev that 



Chandrasekhar 
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must be specified; this may be complicated in the general case. More important 

historically, the equations so developed do not appear to separate. As we'll see 

in a moment, a different approach does yield separable equations, which were 

preferred for much of the history of this field. 

Rather than expanding the metric of the black hole, 

perturbations of its curvature: 



Teukolskvl (|l973l ) examined 



R. 



BH 



+ 5R 



(103) 



The curvature tensor is invariant to first-order gauge transformations, an at- 
tractive feature. This tensor also obeys the nonlinear wave equation (j94p . By 
expanding that equation to linear order in SRa^p,^, Teukolsky showed that per- 
turbations to Kerr black holes are governed by the equation 



2\2 



2 ■ 2 

a sm I 



a?^ - 4 



+- 



4iMamr ^ ^ i ^ 1 



r + ia cos ( 



M(r2 - a2) 



sml 



A 



sm 



m ^ + Aim 



a{r — M) i cos ( 



A 



sm 



^' - (4cot^6' + 2) ^' =T 



(104) 



^ is a complex quantity built from a combinatio i i of c omponents of SRafi/s. 



and describes spacetime's radiation; see iTeukolskyI (119731 ) for details. (We have 



assumed ^ oc e*™"^.) Likewise, T describes a source function built from the 
stress-energy tensor describing a small body orbiting the black hole. Somewhat 
amazingly, Eq. (jl04p separates: putting 



Im 



iTn<f>—iujt 



(105) 



( 19831 ) in fact develops equations to describe Kerr metric perturbations, noting that they rather 



complicated and unwieldy, and are rarely used. 



38 Scott A. Hughes 

and applying a similar decomposition to the source T, we find that Sim{0) is 

a "spin-weighted spheroidal harmonic" (a basis for tensor functions in a non- 
spherical background), and that Rim(r) is governed by a simple ordinary differ- 
ential equation. ^ characterizes Kerr perturbations in much the same way that 
Q [cf. Eq. (j98p ] characterizes them for Schwarzschild. It's worth noting that, 
although the per turbation equations are often so l ved numeric ally, analytic solu- 



tions are known (jMano. Suzuki &: Takasugi 



1993, 



Fiziev 



20091 ). and can be used 



to dramatically improve one's scheme for solving for black hole perturbations. 

Whether one uses this separation or solves Eq. (jl04p directly, solving for per- 
turbations of black holes is now a well-understood enterprise. We now discuss 
how one uses these solutions to model compact binaries. 



4.2 Binary evolution in perturbation theory 

How do we describe the motion of a small body about a black hole? The most 
rigorous approach is to enforce V^T^^, = 0, where T^y describes the small body in 
the spacetime of the large black hole. If we neglect the small body's perturbation 
to the spacetime, this exercise produces the geodesic equation n^V^u'^ = 0, 
where u^^ is the small body's 4- velo city. Geodesic black hole orbits have been 



studied extensively; see, for example. 



Misner. Thorne fc Wheelerl (jl973|), Chapter 



33. A key feature of these orbits is that they are characterized (up to initial 
conditions) by three conserved constants: energy E, axial angular momentum L^, 
and "Carter's constant" Q. If the black hole does not rotate. Carter's constant 
is related to the orbit's total angular momentum: Q{a = 0) = L • L — L^. When 
the black hole rotates rapidly, Q is not so easy to interpret, but the idea that it 
is essentially the rest of the orbit's angular momentum can be useful. 
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Now take into account perturbations from the small body. Enforcing V^T^jy = 

0, we find that the small body follows a "forced" geodesic, 

u'^V^u'' = r , (106) 

where is the covariant derivative in the background black-hole spacetime. The 
novel feature of Eq. (jl06p is the self force f^, a correction to the motion of order 
the small body's spacetime perturbation. The self force is so named because it 
arises from the body's interaction with its own s pacetime correction. 



Self forces have a long pedigree. iDirad (jl938l ) showed that a self force arises 



from the backreaction of an electromagnetic charge on itself, and causes radiative 
damping. Computing the gravitational self force near a black hole is an active 
area of current research. It is useful to break the self force into a dissipative piece, 
/diss' which is asymmetric under time reversal, and a conservative piece, /^ons; 
which is symmetric. These contributions have very different impact on the orbit. 
Dissipation causes t he "conserved" quanti ties {E, L^, Q) to decay, driving inspiral 



of the small body. 



Quinn fc Wald 



(|l999l ) have shown that the rate at which E 



and Lz change due to /^^j^^ is identical to what is found when one computes the 
fluxes of energy and angular momentum encoded by the Isaacson tensor (|72p . 

The conservative self force, by contrast, does not cause orbit decay. "Con- 
served" constants are still conserved when we include this force; but, the orbits 
are different from background geodesies. This reflects the fact that, even neglect- 
ing dissipation, the small body's motion is determined by the full spacetime, not 
just the background black hole. When conservative effects are taken into account, 
one finds that the orbital frequencies are shifted by an amount 

~ X in/M) (107) 
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[where x £ {(j),6,r)]. Because the GWs have spectral support at harmonics of 

the orbital frequencies, these small but non- negligible frequency shifts are directly 

encoded in the waves that the binary gene r ates. Good discussion and a toy model 



can be found in 



Pound. Poisson k Nickell pOOSl l 



To date, not a large amount of work is published regarding self forces and con- 
servative effects for Kerr orbits. There has, ho wever, been enormous progress for 
the case of orbits around non-rotating holes. iBarack &: Sagol (|2007l ) have com- 



pleted an analysis of the full self force for circular orbits about a Schwarzschild 
black hole; generalization to eccentric orbits is in progres s (L. Barack, pr ivate 



communication) . An independent approach developed by iDetweilei 



been found to agree with Barack and Sago extremely well; see 



20081) has 



Sago. Barack &: Detweiler 



(j2008l ) for detailed discussion of this comparison. 



4.3 Gravitational waves from extreme mass ratio binaries 

In this section, we discuss the properties of GWs and GW sources as calculated 
using perturbation theory. As discussed in Sec. ()1.3p . these waves most naturally 
describe extreme mass ratio capture sources. There is also an important overlap 
with the pN results discussed in Sec. EJ By specializing to circular, equatorial 
orbits and considering the limit r S> GM/c^, results from perturbation theory 
agree with pN results for fx/M <^ 1. 

Our goal here is to highlight features of the generic Kerr inspiral waveform. 
As such, we will neglect the conservative self force, which is not yet understood 
for the Kerr case well enough to be applied to these waves. When conservative 
effects are neglected, the binary can be regarded as evolving through a sequence 
of geodesies, with the sequence determined by the rates at which GWs change the 
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"constants" E, L^, and Q. Modeling compact binaries in this limit takes three 

ingredients: First, a description of black hole orbits; second, an algorithm to 

compute GWs from the orbits, and to infer how the waves' backreaction evolves 

us from orbit to orbit; and third, a method to integrate along the or bital sequence 

to bui ld the full waveform. A description of this method is given in 



Hughes et al 



(j2005l ): we summarize the main results of these three ingredients here. 
4.3.1 Black hole orbits. Motion in the vicinity of a black hole can be 
conveniently written in the Boyer-Lindquist coordinates of Eq. (jlOip as r{t), 6{t), 
and (j){t). Because t corresponds to time far from the black hole, this gives a useful 
description of the motion as measured by distant observers. Bound black hole 
orbits are confined to a rGgion near tliG IioIg. They have Tj^j^ 

< r{t) < 

^max and 

^min < d{t) < IT — Omin- Bound orbits thus occupy a torus in th e 3-space near the 



hole's event horizon; an example is shown in Fig. [H taken from 



Drasco Sz Hughes 



(|2006l ). Selecting the orbital constants E, Lz, and Q fully determines r^i^/^^^^^ 
and ^min- It is useful for some discussions to reparameterize the radial motion, 
defining an eccentricity e and a semi-latus rectum p via 

''min — Z ; ) ''max — Z • (^08) 

1 + e 1 — e 



For man y bound black hole orbit s, r{t), 0{t), and (p{t) are periodic (jSchmidt 



see also 



Drasco fc Hughes 



2003; 



20041 ) . (Exceptions are orbits which plunge into the 
black hole; we discuss these below.) Near the hole, the time to cover the full 
range of r becomes distinct from the time to cover the 9 range, which becomes 
distinct from the time to cover 2tt radians of azimuth. One can say that spacetime 
curvature splits the Keplerian orbital frequency O into 0^, and 0^. Figure [2] 
shows these three frequencies, plotted as functions of semi-major axis A for fixed 
values of e and ^min- Notice that all three approach Q oc for large A. 
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4.3.2 Gravitational radiation from orbits. Because their orbits are 

periodic, GWs from a body orbiting a black hole will have support at harmonics 

of the orbital frequencies. One can write the two polarizations 



where 



(109) 



(110) 



The amplitude Hmkn can be found by solving the Teukolsky equation (|104p 
using the decompositio n (jlOSp : details for the general case can be found in 



Drasco &: Hughes! (j2006l ). An example of a wave from a geodesic orbit is shown 
in Fig. [3l Note the different timescales apparent in this wave; they are due to the 
three distinct frequencies of the underlying geodesic orbit (and their harmonics). 

The expansion (|109p does not work well for orbits that plunge into the black 
hole; those orbits are not periodic, and cannot be expanded using a set of real 
frequencies. A better way to calculate those waves is to solve the Teukolsky 
equation (|104p without introducing the decomposion (jlOSp . Results for waves 
from plunging orbits in the l angu a ge of perturbat i on th eory was first given by 



Nagar. Damour fc Tartaglia 



(120071) 



SundararaianI (|2008l ) has recently extended 



the cases that we can model to full generality. 

As mentioned in Sec. 14. 2^ it is fairly simple to compute the flux of energy E 
and angular momentum from the Isaac son tensor, Eq. (j72p . once the waves 



are known. Recent work (^Ganz et al 



20071 ) has shown that a similar result de- 



scribes Q. Once E, L^, and Q are known, it is straightforward to evolve the 
orbital elements r-jnin/^nax and 6'min, specifying the sequence of orbits through 
which gravitational radiation drives the system. Figure U] gives an example of 
how orbits evolve when their eccentricity is zero. 
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4.3.3 Evolving THROUGH AN ORBITAL SEQUENCE. It is not too difficult to 

compute the sequence of orbits [parameterized as E{t), Lz{t), Q{t) or ^jnin/jnaxl*)) 

Ommif)] that an inspiraUng body passes through before finally plunging into its 

companion black hole. Once these are known, it is straightforward to build the 

worldline that a small body follows as it spirals into the black hole. From the 

worldline, we can build a source function T{t) for Eq. (jl04p and compute the 

evolving inspiral waves. Figure [5] gives an exa mple of a wave arisin g from , the 



Sundararajan et al 



((20081) for 



inspiral of a small body into a black hole; see 
details of how these waves are computed. 

5 Post-Newtonian theory 

Suppose we cannot use mass ratio as an expansion parameter. For instance, if 
the members of the binary are of equal mass, then rj = mim2/{mi +771,2)^ = 0.25. 
This is large enough that neglect of 0{rf') and higher terms is problematic. The 
techniques discussed in Sec. H] will not be appropriate. 

If the mass ratio is not a good expansion parameter, the potential (f) = GMjrf? 
may be. The post-Newtonian (pN) expansion of GR results when we use (j) as our 
expansion parameter. We now summarize the main concepts which underlie the 
pN formalism, turning next to a discussion of the pN w a veform and its interesting 



features. Much of our discussion is based on 



Blanched ((20061). 



5.1 Basic concepts and overview of formalism 

One typically begins the pN expansion by ex amining the Ein stein field equations 
in harmonic or deDonder coordinates (e.g., iWeinberei Il972l . Sec. 7.4). In these 
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coordinates, one defines 

h^'^V^g^'-v'"' , (111) 

where g is the determinant of Qf^u- This looks similar to the flat spacetime per- 
turbation defined in Sec. 13.11 however, we do not assume that h is small. We 
next impose the gauge condition 

= . (112) 

With these definitions, the exact Einstein field equations are 

□/i"^ = l^r°^ , (113) 

where □ = rj^'^dadfs is the flat spacetime wave operator. The form of Eq. (|113p 
means that the radiative Green's function we used to derive Eq. ()4ip can be 
applied here as well; the solution is simply 

c J |x — x'l 

Formally, Eq. (|114p is exact. We have swept some crucial details under the 
rug, however. In particular, we never defined the source r"^. It is given by 

= {-g)T^P + —— , (115) 

where T"^ is the usual stress energy tensor, and A"^^ encodes much of the non- 
linear structure of the Einstein field equations: 

A"^ = lQ^l{-g)tll + ^^h'^^'^f,h'^'' -df.d^h'^'^h''^ (116) 
= A^^'^f/i,/!] +M"^[/i,/i,/i] +L"^[/i,/i,/i,/i] +0(^1^) . (117) 

On the first line, t'^ is the Landau-Lifshitz pseudotensor, a quantity which (in 
certain gauges) allows us to describe how GWs carry energy through spacetime 



( Landau &: Lifshitz 



19751 . Sec. 96). On the second line, the term N°^^[h, h] means 
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a collection of terms quadratic in h and its derivatives, M"^[h, h, h] is a cubic 

term, etc. Our solution h"^ appears on both the left- and right-hand sides of Eq. 

(jll4p . Such a structure can be handled very well iteratively. We write 

oo 

h'^P = Y^G''hf . (118) 

n=l 

The n = 1 term is essentially the linearized solution from Sec. 13.11 To go higher, 
let An^ denote the contribution of A"^ to the solution hn^ . We find 

Af = N''('[hi,hi], (119) 

Af = M"^[/ii, /ii, /ii] + iV"^[/i2, ^i] + iV"^[/ii, /i2] , (120) 

etc.; higher contributions to A"'' can be found by expanding its definition and 
gathering terms. By solving the equations which result from this procedure, it 
becomes possible to build the spacetime metric and describe the motion of the 
members of a binary and the radiation that they emit. 

5.2 Features of the post-Newtonian binary waveform 

The features of the pN binary waveform are most naturally understood by first 
considering how we describe the motion of the members of the binary. Take those 
members to have masses mi and m2, let their separation be r, and let f point 
to body 1 from body 2. Then, in the harmonic gauge used for pN theory, the 
acceleration of body 1 due to the gravity of body 2 is 

a. = ao + ^2 + + + + a.7 . . . . (121) 

The zeroth term, 

ao = --^r, (122) 
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is just the usual Newtonian gravitational acceleration. Each a„ is a pN correction 

of order {v/cY\ The first such correction is 



a2 



5G^mim2 AG'^ml ^ Gm2 f 3 



+ 



+ 



^(f • va)^ -?;? + 4vi • V2 -2vi 



r 

72 



(123) 



For the acceleration of body 2 due to body 1, exchange labels 1 and 2 and replace 
f with — f . Note that a2 changes the dependence of the acceleration with respect 
to orbital separation. It also shows that the acceleration of body 1 depends on 
its mass mi. This is a pN manifestation of the "self force" discussed in Sec. 14.21 
So far, the pN acceleration has been computed to order {v/cY. As we go to 
high order, the expre s sions for a„ become quite lengthy. An excellent summary 



is given in iBlanchetl (|2006l ). Eq. (131) and surrounding text. (Note that the 
expression for ag fills over two pages in that paper!) 

At higher order, we also find a distinctly non-Newtonian element to binary 
dynamics: its members spins precess due to thei r motion in the binary' s curved 



spacetime. If the spins are Si and S2, one finds (jThorne &: Hartle 



19851 ^ 



dSi 



G 
G 



dt c^r'^ 



3 m2 , 
2 + 77— /xVMrL 

2 mi 

2 m2 



xSi + 



XS2 + 



G 



G 



is2-^(S2-L)L 



iSi-^(Si.L)L 



xSi , (124) 
XS2 . (125) 



We now discuss the ways in which aspects of pN binary dynamics color a 
system's waves. 

5.2.1 Gravitational- WAVE amplitudes. Although a binary's dominant 
waves come from variations in its mass quadrupole moment, Eq. (|54p shows us 
that higher moments also generate GWs. In the pN framework, these moments 
contribute to the amplitude of a binary's waves beyond the quadrupole form, Eq. 
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(j89p . Write the gravitational waveform from a source as 

2GM /irGMfY^^ 
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(126) 



where v = {ttGM f /c^y^^ is roughly the orbital speed of the binary's members 
(normalized to c). The contributions ^ reproduce the waveform presentee 
i n Eg . (j89p . The higher-order terms and i^l x can be found in iBlanchetl 

(|2006l ). his Eqs. (237) through (241). For our purposes, the key point to note is 
that these higher-order terms introduce new dependences on the binary's orbital 
inclination and its masses. As such, measurement of these terms can provide 
additional constraints to help us understand the system's characteristics. Figure 
[6] illustrates the three contributions Hq, H1/2, and Hi to a binary's GWs. 
5.2.2 Orbital phase. The motion of a binary's members about each other 
determines the orbital phase. Specializing to circular orbits, we can determine 
the orbital frequency from the acceleration of the binary's members; integrating 
up this frequency, we define the binary's phase ^{t). The first few terms of this 



phase are given by ^Blanchet et al 



19951 ^ 



5GM 



1 5/8 



1 + 



/3715 , 55 n 
V8064 



/ 9275495 284875 u 1855 15 , 

+ \ — H At H crit) Q~^'^ 

V 14450688 258048 M 2048 64 ^7 



where 



G 



c r] 



{tc - t) . 



(127) 



(128) 



5GM 

Notice that the leading term is just the Newtonian quadrupole phase, Eq. ([90]) . 
Each power of connects to a higher order in the pN; Eq. (jl27p is taken to "sec- 
ond post-Newtonian" order, which means that corrections of {vl c)^ are included. 



Corrections to order (v/c)^ are summarized in iBlanchetl ( 20061 ). In addition to 



the chirp mass M, the reduced mass fi enters $ when higher order terms are 
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included. The high pN terms encode additional information about the binary's 

masses. At least in principle, including higher pN effects in our wave model makes 

it possible to determine both chirp mass and reduced mass, fully constraining the 

binary's masses. 

Equation (jl27p also depends on two parameters, (3 and a, which come from 
the binary's spins and orbit orientation. The "spin-orbit" parameter /3 is 



1 



i=l 



113 (^) +75, 



2 ' 

mf 



(129) 



the "spin-spin" parameter a is 



V 



(Blanchet et al 



721(L-Si)(L-S2) - 247Si • Sa 



(130) 



19951 ). As we'll see in Sec. [H these parameters encode valuable 



information, especially when spin precession is taken into account. 

The <^ M limit of Eq. ()127p can be computed with black hole perturbation 
theory (Sec. H]) evaluated for circular orbits with r ^ GM/c^. The orbital phase 
is found by integrating the orbital frequency. By changing variables, one can 
relate this to the orbital energy and the rate at which GWs evolve this energy: 

poih 



(131) 



The orbital energy E""^^ is simple to calculate and to express as a function of 
orbital frequency. For example, for orbits of non-rotating black holes, we have 

1 - 2v^/c^ 



v/r^3Z7cLl. 



(132) 



where v = rO,. For circular, equatorial orbits. iMino et al 



(|l997l ) analytically solve 



the Teukolsky equation as an expansion in v, calculating dE^^ /dt to 0[{v/cf^]. 
This body of work confirms, in a completely independent way, all of the terms 
which do not depend on the mass ratio /i/M in Eq. (jl27p . The fact that these 
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terms are known to such high order is an important input to the effective one- 
body approach described in Sec. 15.31 

5.2.3 Spin precession. Although the spin vectors Si and S2 wiggle around 
according to the prescription of Eqs. ()124p and (jl25p . the system must preserve 
a notion of global angular momentum. Neglecting for a moment the secular 
evolution of the binary's orbit due to GW emission, pN encodes the notion that 
the total angular momentum 

J = L + Si + S2 (133) 

must be conserved. This means L must oscillate to compensate for the spins' 
dynamics, and guarantees that, when spin precession is accounted for in our evo- 
lutionary models, the phase parameters P and a become time varying. Likewise, 
the inclination angle l varies with time. Precession thus leads to phase and am- 
plitude modulation of a source's GWs. Figure [7] illustrates precession's impact, 
showing the late inspiral waves for binaries that are identical aside from spin. 

5.3 The effective one-body approach 

Because pN techniques are based on an expansion in = GM/rc^, it had been 
thought that they would only apply for r ^ lOGM/c^ , and that numerical rela- 
tivity would be needed to cover the inspiral past that ra dius, through the fina. 



plunge and merger. This thinking was radically changed bv lBuonanno &: Damour 



(|l999l ). which introduced the effective one-body approach to two-body dynamics. 
This technique has proved to be an excellent tool for describing the late inspi- 
ral, plunge, and merger of two black holes. We now describe the key ideas of 
this app r oach; our description owes much to the helpful introductory lectures by 



Damoui] ^2008) . 
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As the name suggests, the key observation of this approach is that the motion 

of two bodies (mi, 777-2) about one another can be regarding as the motion of a 
single test body of mass /i = mim2/{mi + 1712) in some spacetime. One begins 
by examining the Hamiltonian which gives the conservative contribution to the 
equations of motion. Let the binary's momenta be pi^2 and its generahzed po- 
sitions qi,2- If we work in the center of mass frame, then the Hamiltonian can 
only be a function of the relative position, q = qi — q2, and can only depend 
on the momentum p = Pi = — P2- For example, the conservative motion can be 
described to second-post-Newtonian order [i.e., 0{v^/c^)] with the Hamiltonian 

H{ci, p) = Ho{p, q) + \h2{p, q) + ^H^ip, q) , (134) 

where Hq = |pp/2/x + GMfi/\q\ encodes the Newtonian dynamics, and ^^2,4 de- 
scribes pN corrections to that motion. A binary's energy and angular momentum 
can be found from this Hamiltonian without too much difficulty. 
The next step is to write down an effective one-body metric, 

ds'^ = -A{R)c^dT'^ + B{R)dR^ + R^idO"^ + sin^ 9 dcf?) , (135) 

where A{R) = 1 + ai{GM / Rc^) + a2{GM / Rc^Y + • • •! a similar expansion de- 
scribes B{R). The coefficients Oi depend on reduced mass ratio, r] = fi/M. The 
effective problem is then to describe the motion of a test body in the space- 
time (|135p . By asserting a correspondance between certain action variables in 
the pN framework and in the effective framework, the coefficients Oi are com- 
pletely fixed. For example, one finds that, as r/ — > 0, the metric (jl35p is simply 
the Schwarzschild spacetime. The effective problem can thus be regarded as the 
motion of a test body around a "deformed" black hole, with rj controlling the 



deformation. See 



Damouij (|2008l ) and references therein for further discussion. 



Compact binary gravitational waves 51 
One must also describe radiation reaction in the effecti ve one-body approach. A 



key in novation introduced by Damour and coworkers [see 



Damour. Iyer &: Sathvaprakash 



(|l998l ) for the original presentation of this idea] is to re-sum the pN results for 
energy loss due to GWs in order to obtain a result that is good into the strong 
field. In more detail, we put 

The function J^a is known to rather high order in orbital velo city v hy a combi- 



nation of analyses in both pN theory [see, e.g., 
to analytic expansion of results in perturbation theory 
be written 

:Fiv) - ^2^-^4o5 



Blanchetl (120061) for a rev iew] and 



Mino et al 



19971 ). It can 



where 



5c5 



F(z;) = l + ^a„Q"^' + ^6„log(7;/c)(^)'' . (138) 
Post-Newtonian theory allows us to compute ctn including co ntributions in n/M, 



(137) 



n/2 



up to n = 7, a nd shows that bn / for n = 6 



Blanchet 



bation theory 



Mino et al 



20061 . Eq. (168)]. Pertur- 



Il997l . Eq. (4.18)] gives us the 0[{n/Mf] contributions 
for ttn up to n = 11, and shows that bn ^ for n = 8, 9, 10, 11. 

The resummation introduced by Damour, Iyer & Sathyaprakash requires fac- 
toring out a pole at V = V in the function F{v) and then reorganizing the remain- 
ing terms using a Fade approximant: 



F'^v) = (1 - v/vy^ P [(1 - v/v)F{v)] 



(139) 



The approximant P converts an A^-th order polynomial into a ratio of A/2-th 
order polynomials whose small v expansion reproduces the original polynomial: 



P 



N 



1 + y^Cn(t;/c) 



n=l 



l + EnLleniv/cY 



(140) 
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Using this approach to define the evolution of a system due to GW backreaction, 

it is not so difficult to compute the waves that a binary generates as its members 

spiral together. Indeed, by augmenting these waves with the "ringdown" that 

comes once the spacetime is well-described by a single black hole, the effective 

one-body approach has recently had great success in matching to the waveforms 

that are produced by numerical relativity simulations. We defer a discussion of 

this matching until after we have described numerical relativity in more detail, 

in order that the effectiveness of this comparison can be made more clear. 

6 Numerical relativity 

Numerical relativity means the direct numerical integration of the Einstein field 
equations, evolving from an "initial" spacetime to a final state. This requires re- 
thinking some of our ideas about GR. As a prelude, consider Maxwell's equations, 
written in somewhat non-standard form: 

V • E = 47rp , V • B = ; (141) 

— = -cV X E , — = 47rJ - cV X B . (142) 

dt dt ^ ^ 

These equations tell us how E and B are related throughout spacetime. Notice 
that Eqs. ()14ip and ()142p play very different roles here. The divergence equations 
contain no time derivatives; if we imagine "slicing" spacetime into a stack of 
constant time slices, then Eq. (|14ip tells us how E and B are constrained on 
each slice. By constrast, the curl equations do include time operators, and so 
tell us how E and B are related as we evolve from slice to slice. We turn now to 
developing the Einstein equations into a form appropriate for evolving from an 
initial form; our discussion has been heavily influenced by the nicely pedagogical 
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presentation of iBaumgarte &: Shapiro 



(120031). 
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6.1 Overview: From geometric equations to evolution equations 

How do we use Eq. (|29p to evolve spacetime from some initial state? The Einstein 
field equations normally treat space and time democratically — no explicit notion 
of time is built into Eq. (j29p . The very question requires us to change our thinking: 
"evolving" from an "initial" state requires some notion of time. 

Suppose that we have chosen a time coordinate, defining a way to slice space- 
time into space and time. We must reformulate the Einstein field equations using 
quantities defined solely on a given time slice. Figure [5] illustrates how two nearby 
time slices may be embedded in spacetime. Once time is set, we can freely choose 
spatial coordinates in each slice; we show on both slices. 

Let n be normal to the bottom slice. The lapse a = dr/dt sets the proper time 
experienced by an observer who moves along n; the shift /?* tells us by how much 
X* is displaced ("shifted") on the second slice relative to the normal observer. We 
will soon see that a and /?* are completely unconstrained by Einstein's equations. 
They let us set coordinates as conveniently as possible, generalizing the gauge 
generator used in linearized theory (cf. Sec. [3|) to the strong field. 

The proper spacetime separation of and + dx* is then 

ds2 ^ _a'^dt^ + g..(^dx' + P'dt){dx^ + P^dt) . (143) 

(In this section, we will put c = 1; various factors become rather unwieldy other- 
wise.) We now have a form for the metric of spacetime, a notion of constant time 
slices, and the normal to a slice n. Because we are interested in understanding 
quantities which "live" in a given slice (i.e., orthogonal to the normal n), we build 
the projection tensor j^^, = g^^ + n^rii,. This tensor is just the metric for the 
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geometry in each slice. We can choose coordinates so that ^tt = 1u = 0, and 

7ij = dij] we will assume these from now on. 

We now have enough pieces to see how to build the field equations in this 

formalism: We take Eq. (j29p and project components parallel and orthogonal to 

n. Consider first the component that is completely parallel to n: 

Gapn'^n^ = SirGTapn'^nf^ — > R + K^ - KijK'^ = WnGp . (144) 



[See iBaumgarte Sz Shapiro! (120031 ) for a detailed derivation of Eq. ()144p .] In this 



equation, R is the Ricci scalar for the 3-metric jij, p = Tajsn'^n^, and 

Kij = -'Yi'^j/Vanp 

= ^{-daij + D,f)j+Djl3i) (145) 

is the extrinsic curvature. (The operator Di is a covariant derivative for the 
metric "fij-) It describes the portion of the curvature which is due to the way 



that each constant time slice is embedded in t he full spacetime. Equation )T1^ 
is known as the Hamiltonian constraint. [See iBaumgarte fc Shapirol (|2003l ) for 



details of how to go from the first line of (jl44p . which is a definition, to the second 
line, which is more useful here.] Notice that it contains no time derivatives of 
Kij. This equation is thus a constraint, relating data on a given timeslice. 
Next, components parallel to n on one index and orthogonal on the other: 

G„/3n%^ = 87rGr„^n%^ D^K^ - DiK = SnGji (146) 

The matter current ji = —T^pn'^^i^ ■ Equation (|146p is the momentum constraint] 
notice it also has no time derivatives of Kij. 
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Finally, project completely orthogonal to n: 

-DiDjO + a Rij - 2KikK^j + KKij - SvrGa (matter) 
(i^DkKij + K,kDj(3'' + KkjD,f3'' . (147) 



dtKij 



[We have abbreviated a combination of projections of t he stress-energy tensor 



as "matter." Interested readers can find more details in 



Baumgarte &: Shapiro 



(|2003l ).] Combining Eqs. (jl47p with (jl45p gives us a full set of evolution equations 
for the metric and the extrinsic curvature, describing how the geometry changes 
as we evolve from time slice to time slice. 



The field equations sketched here are the A DM equations (jArnowitt. Deser &: Misner 



19621 ) . Today, most groups work with modified versions of these equations; a par 



(|l999l ) , building on foundational work by 



em, developed by 


1„ — „ 1 — * 

Baumgarte &: ShaDiro 


Shibata & Nakamura 


(1995 


). InBSSN, 



one rewrites the spacetime metric as 



(148) 



where (p is chosen so that e^^*^ = det(7jj). With this choice, det(7jj) = 1. Roughly 
speaking, the decomposition (jl48p splits the geometry into "transverse" and "lon- 
gitudinal" degrees of freedom (encapsulated by 7jj and 4>, respectively.) One sim- 
ilarly splits the extrinsic curvature into "longitudinal" and "transverse" parts by 
separately treating its trace and its trace-free parts: 

1 



Aij — Kij gTij-^ • 



It is convenient to conformally rescale Aij, using A. 



(149) 



e-'^'I'AiA. One the n 



develo ps evolution equations for cj), jij, K, and Aij. See iBaumgarte &: Shapiro 



(jl999l ) for detailed discussion. 
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6.2 The struggles and the breakthrough 

Having recast the equations of GR into a form that lets us evolve forward in 
time, we might hope that simulating the merger of two compact bodies would 
now not be too difficult. In addition to the equations discussed above, we need a 
few additional pieces: 

1. Initial data; i.e., a description of the metric and extrinsic curvature of a 
binary at the first moment in a simulation. Ideally, we might hope that 



rated bodies, (e.g.. 



Nissanke 


2006 


, and 


Yunes et al. 


2006) 



20061 ). However, any 



method which can produce a bound binary with specified masses mi, 772-2 
and spins Si, S2 should allow us to simulate a binary (although it may be 
"contaminated" by having the wrong GW content at early times). 

2. Gauge or coordinate conditions; i.e., an algorithm by which the lapse a 
and shift /3* are selected. Because these functions are not determined by 
the Einstein field equations but are instead freely specified, they can be 
selected in a way that is as convenient as possible. Such wonderful freedom 
can also be horrible freedom, as one could choose gauge conditions which 
obscure the physics we wish to study, or not facilitate a stable simulation. 

3. Boundary conditions. If the simulation contains black holes, then they 
should have an event horizon from which nothing comes out. Unfortu- 
nately, we cannot know where horizons are located until the full spacetime 
is built (though we have a good estimate, the "apparent horizon," which can 
be computed from information on a single time slice). They also contain 
singularities. Hopefully, event horizons will prevent singular fields from 
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contaminating the spacetime. The computation will also have an outer 

boundary. Far from the binary, the spacetime should asymptote to a flat 

(or Robertson- Walker) form, with a gentle admixture of outgoing GWs. 

How to choose these ingredients has been active research for many years. Early 
on, some workers were confident it was just a matter of choosing the right com- 



bination of methods and the pieces would fall into place. 
was nicely encapsulated by the following quote from 



he ini tial optimism 



Thornd (119871 1 . p. 379: 



. . . numerical relativity is likely to give us, in the next five years or 
so, a detailed and highly reliable picture of the final coalescence and 
the wave forms it produces, including the dependence on the holes' 
masses and angular momenta. 

For many years, Thome's optimism seemed misplaced. Despite having a de- 
tailed understanding of the principles involved, it seemed simply not possible to 
evolve binary systems for any interesting length of time. In most cases, the bi- 
nary's members could complete a fraction of an orbit and then the code would 
crash. As Joan Centrella has emphasized in several venues, by roughly 2004 "Peo- 
ple were beginning to say 'numerical relativity cannot be done'" (J. Centrella, 
private communication) . 

A major issue with many of these simulations appeared to be constraint violat- 
ing modes. These are solutions of the system {dtjij, dtKij) that do not satisfy the 
constraint equations (jl44p and (jl46p . As with the Maxwell equations, one can 
prove that a solution which satisfies the constraints initially w ill continu e to sa t- 
isfy them at later times in the continuum limit [cf. Sec. IIIC of [ 
Unfortunately, numerical relativity does not work in the continuum limit; and, 
even more unfortunately, constraint violating modes generically tend to be un- 



Pretoriusl ([20091)]. 
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stable (|Lindblom &: Scheel 
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20021 ). This means that small numerical errors can 
"seed" the modes, which then grow unstably and swamp the true solution. The 
challenge was to keep the seed ins tabilities as small as possible and th en to prevent 



them from growing. In this way, 



Briigmann. Tichv &: JansenI (|2004l ) were able to 



compute one full orbit of a black hole bi nary before thei r simulation crashed. 



It was thus something of a shock when 



Pretoriud (|2005l ) demonstrated a binary 



black hole simulation that executed a full orbit, followed by merger and ringdown 
to a Kerr black hole, with no crash apparent. Pretorius used a formulation of the 
Einstein equations based on coordinates similar to the de Bonder coordinates 
des cribed in Sec. [5l These are known as "generalized harmonic coordinates"; 



see iPretoriusI (j2009l ) for detailed discussion. Because his success came with such 
a radically different system of equations, it was suspected that the ADM-like 
equations might have to be abandoned. These concerns were allayed by near 
simultaneous discoveries by numerical relativity groups then at the University 



of Texa s in Browsnville (iCampanelli et al. 



Center (Baker et al 



20061 ^ and the Goddard Space Fhght 



20061 ). The Campanelli et al. and Baker et al. groups both 
used the BSSN formalism described at the end of Sec. 16.11 To make this approach 
work, they independently discovered a gauge condition that allows the black holes 
in their simulations to move across their computational grid. (Earlier calculations 
typically fixed the coordinate locations of the black holes, which means that the 
coordinates effectively co-rotated with the motion of the binary.) It was quickly 
shown that this approac h yielded results that agreed excellently with Pretorius' 



setup (IBaker et al. 



2003). 

Implementing the so-called "moving puncture" approach was sufficiently sim- 
ple that the vast majority of numerical relativity groups were able to immedi- 
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ately apply it to their codes and begin evolving binary black holes. These tech- 



2008 



(Shibata & Urvu 


2006 


Etienne et al. 


2008 


Baiotti. Giacomazzo & Rezzolla 



Shibata et al 



20091 ). In the past four years, we have thus moved from a 
state of being barely able to model a single orbit of a compact binary system in 
full GR to being able to model nearly arbitrary binary configurations. Though 
Thome's 1987 prediction quoted above was far too optimistic on the timescale, 
his prediction for how well the physics of binary coalescence would be understood 
appears to be exactly correct. 

6.3 GWs from numerical relativity and effective one body 

Prior to this breakthrough, the effective one-body approach gave the only strong- 
field description of GWs from the coalescence of two black holes. Indeed, these 
techniques made a rather strong prediction: The coalescence waveform should be 
fairly "boring," in the sense that we expect the frequency and amplitude to chirp 
up to the point at which the physical system is well modeled as a single deformed 
black hole. Then, it should rapidly ring down to a quiescent Kerr state. 

Such a waveform is indeed exactly what numerical simulations find, at least 
for the cases that have been studied so far. It has since been found that pre- 
dictions from the effective one-body formalism give an outstanding description 
of the results from numerical relativity. There is some freedom to adjust how 
one matches effective one-bod y waves to the nu i nerica l rela tivity output [e . g., the 



choice of pole v in Eq. (jl39p ]: 



Buonanno et al 



(120071) and 



D amour et al 



( 2008 



describe how to do this matching. 

Figure m taken from Buonanno et al., gives an example of how well the wave- 
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forms match one another. Over the entire span computed, the two waveforms 

differ in phase by only a few hundredths of a cycle. The agreement is so good that 

one can realistically imagine "calibrating" the effective one-body waveforms with 

a relatively small number of expensive numerical relativity computations, and 

then densely sampling the binary parameter space using the effective one-body 

approach. 

7 Gravitational-wave recoil 

That GWs carry energy and angular momentum from a binary, driving its mem- 
bers to spiral together as described in Sec. 13. 3^ is widely appreciated. Until 
recently, it was not so well appreciated that the waves can carry linear momen- 
tum as well. If the binary and its radiation pattern are asymmetric, then that 
radiation carries a net flux of momentum given by 

^ = ^ [ dnr'^'n', (150) 



dt 

where T^^ is the energy- flux component of the Isaacson tensor (I72p . n* is the i-th. 
component of the radial unit vector, and the integral is taken over a large sphere 
(R — > oo) around the source. Recent work has shown that the contribution 
to this "kick" from the final plunge and merger of coalescing black holes can 
be particularly strong. We now summarize the basic physics of this effect, and 
s urvey recent result s. 



BekensteinI (jl973l ) appears to have first appreciated that the momentum flux 



FitchettI (119831 ) then 



of GWs could have interesting astrophysical consequences, 
estimated the impact this flux could have on a binary. An aspect of the problem 
which Fitchett's analysis makes very clear is that the recoil comes from the beat- 
ing of different multipolar contributions to the recoil: If one only looks at the 
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quadrupole part of the GWs [cf. Eq. (j54p ]. the momentum flux is zero. Fitchett's 

analysis included octupole and current-quadrupole radiation, and found 

c 1450 km/sec M(^|W^:!)\ (151) 

/max \ -Ktorm / 

where f{q) = g^(l — q)/{l + gives the dependence on mass ratio q = raxjra^. 
This function has a maximum at g ~ 0.38. The radius i?torm describes when wave 
emission cuts off; for systems containing black holes this will scale with the total 
mass. Thus, the recoil does not depend on total mass, just mass ratio. 

Fitchett's analysis is similar to our discussion in Sec. 13.31 in that Newtonian 
dynamics are supplemented with multipolar wave emission. Because the effect 
is strongest when -Rterm is smallest, it was long clear that a proper relativis- 
tic analysis was needed to get this kick correct. Indeed, a prescient analysis 



by 



Redmount &: Rees 



(|l989l ) suggested that binaries containing rapidly spinning 
black holes were likely to be especially interesting; as we shall see in a few mo- 
m ents, they were absolute l y corr ect . 



Favata. Hughes fc Hold (|2004l ) provided the first estimates of recoil which did 
the strong-field physics more-or-less correctly. They used the perturbative tech- 
niques described in Sec. [H arguing that one can extrapolate from the small mass 
ratio regime '10 q^ 0.2 or so with errors of a few tens of percent. Unfortunately, 
their code at that time did not work well with plunging orbits, so they had large 
error bars. They did find, however, that the maximum recoil probably fell around 
fkick — (250 ± 150) km/se c, at least if no more than one bo dy spins, and if the 



spin and orbit are aligned. iBlanchet. Qusailah &: Willi (|2005l ) revisited this treat 



ment, with a particular eye on the final plunge waves, using the pN methods 
outlined in Sec. [5l Their results were consistent with Favata et al., but reduced 
the uncertainty substantially, finding a maximum recoil ^kick — (220ib50) km/sec. 
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These numbers stood as the state-of-the-art in black hole recoil for several 

years, until numerical relativity's breakthrough (Sec. [6]) made it possible to study 

black hole mergers without any approximations. The Favata et al. and Blanchet 

et al. numbers turn out to agree quite well wit h predictions for the merger of 



non-spinning black holes; 



Gonzalez et al 



(|2007bl ) find the maximum kick for non- 



spinning merger comes when q = 0.38, yielding ?;kick = (175 it 11) km/sec. When 
spin is unimportant, kicks appear to be no larger than a few hundred km/sec. 



When spin is importaii t, the kick can be subs t antia. 



by 



Gonzalez et al 



(|2007al l and 



Campanelli et al 



ly larger. Recent work 



(|2007l ) shows that when the 



holes have large spins and those spins are aligned just right (equal in mag- 
nitude, antiparallel to each other, and orthogonal to the orbital angular mo- 
mentum), the recoil can be a few thousand km/sec. Detailed parameter ex- 
ploration is needed to assess how much of this maximum is actually achieved; 
early work on this problem is finding that large kicks (many hundreds to a few 
thousand of km/se c) can be achieved for va, r ious spin orientatio ns as long as 



th e spins are large (iTichv &: 



by 



Marronett 



Bovle. Kesden &: Nissanke 



20071 . 



Pollnev et al. 



20071 ): recent work 



20081 ) shows how the recoil can depend on spin 
and spin orientation, suggesting a powerful way to organize the calculation to 
see how generic such large kicks actually are. That the maximum is so much 
higher than had been appreciated suggests that substantial recoils may be more 
common than the Favata et al. and Blanchet et al. calculations led us to expect 



(jSchnittman &: Buonanno 



2003). 

Many recent papers have emphasized that kicks could have strong astrophysical 
implications, ranging from escape of the black hole from its host galaxy to shocks 



in material accreting onto the large black hole. The first possible detection of 
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black hole recoil was recently announced (|Komossa. Zhou &: Luii2008i ). As that 

claim is assessed, we anticipate much activity as groups continue to try to identify 

a signature of a recoiling black hole. 

8 Astronomy with gravitational waves 

Direct GW measurement is a major motivator for theorists seeking to understand 
how binary systems generate these waves. The major challenge one faces is that 
GWs are extremely weak; as we derive in detail in this section, a wave's strain 
h sets the change in length AL per length L in the arms of a GW detector. 
Referring to Eq. (j89p . we now estimate typical amplitudes for binary sources: 



h 



2GM f7rGMf\ 



amp - ^2 2) 



2/3 



c3 ) 



^P_23 ^ 2.8 Mo V / ^ V Mpc 



M J viooHzy V D 

ii^i^X'' (fy" f ^ I . (152) 



M J V 10-3 Hz 7 V D 
On the last two lines, we have specialized to a binary whose members each 
have mass M/2, and have inserted fiducial numbers corresponding to targets 
for ground-based detectors (second line) and space-based detectors (third line). 

In the remainder of this section, we summarize the principles behind modern 
interferometric GW detectors. Given that these principles are likely to be novel 
for much of the astron omical community, we present this material in some depth. 
Note that 



Finn! (j2009l ) has recently flagged some important issues in the "stan- 
dard" calculation of an interferometer's response to GWs (which, however, do not 
change the final results). For the sake of brevity, we omit these issues and rec- 
ommend his article for those wishing a deeper analysis. We then briefly describe 
existing and planned detectors, and describe how one measures a binary's signal 
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with these instruments. This last point highhghts why theoretical modeling has 

been so strongly motivated by the development of these instruments. 



8.1 Principles behind interferometric GW antennae 

As a simple limit, treat the spacetime in which our detector lives as flat plus a 
simple GW propagating down our coordinate system's z-axis: 

ds"^ = -c^dt^ + (1 + h)dx^ + (1 - h)dy'^ + dz'^ , (153) 

where h = h{t — z). We neglect the influence of the earth (clearly important 
for terrestrial experiments) and the solar system (which dominates the spacetime 
of space-based detectors). Corrections describing these influences can easily be 
added to Eq. ()153p : we neglect them, as they represent influences that vary on 
much longer timescales than the GWs. 

Figure [10] sketches an interferometer that can measure a GW. Begin by ex- 
amining the geodesies describing the masses at the ends of the arms, and the 
beam splitter at the center. Take these objects to be initially at rest, so that 
((ix^/dr)before = (c, 0,0,0). The GW shifts this velocity by an amount of or- 
der the wave strain: ((ix^/dr)after = {dx^ / dr^eiove + 0{h). Now examine the 
geodesic equation: 

d^x^ ^- dx'^ dx^ 

All components of the connection are 0{h). Combining this with our argument 
for how the GW affects the various velocities, we have 

d'^x^ dx^ dx'^ ,^,,0. , . 

-r^ + r^oo-^-^ + 0(/i2) = 0. (155) 
ar^ dr or 

Now, 

r^oo = ^g^'' {dogw + dogok - dtgoo) = (156) 
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as the relevant metric components are constant. Thus, 
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(157) 



The test masses are unaccelerated to leading order in the GW amplitude h. 

This seems to say that the GW has no impact on the masses. However, the 
geodesic equation describes motion with respect to a specified coordinate system. 
These coordinates are effectively "comoving" with the interferometer's compo- 
nents. This is convenient, as the interferometer's components remain fixed in our 
coordinates. Using this, we can show that the proper length of the arms does 
change. For instance, the x-arm has a proper length 



^ V9^dx = J^ VTThdx^j^ {^^ + dx = L (^1 + . (158) 



Likewise, the y-arm has a proper length Dy = L{1 — h/2). 

This result tells us that the armlengths as measured by a ruler will vary with h. 
One might worry, though, that the ruler will vary, cancelling the measurement. 
This does not happen because rulers are not made of freely-falling particles: The 
elements of the ruler are bound to one another and act against the GW. The ruler 
will feel some effect from the GW, but it will be far smaller than the variation in 
the separation. 

The ruler used by the most sensitive current and planned detectors is based 
on laser interferometry. We now briefly outline how a GW imprints itself on 



the phase o f the int e rferom eter sketched in Fig. TOl 



recommend 



or further discussion, we 



FinnI (j2009l ) for more detailed 



Faraonil (j2007l ): we also recommend 
discussion. In the interferometer shown, laser light enters from the left, hits 
the beam splitter, travels down both the x- and y-arms, bounces back, and is 



recombined at the beam splitter. Begin with light in the x-arm, and compute 
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the phase difference between hght that has completed a round trip and light that 

is just entering. The phase of a wavefront is constant as it follows a ray through 

spacetime, so we write this difference as 

A^x = $(rround-trip) " ^(0) = ^proper X (proper round-trip travel time) , (159) 

where Wproper is the laser frequency measured by an observer at the beam splitter. 
Detecting a GW is essentially precision timing, with the laser acting as our clock. 

Consider first the proper frequency. Light energy as measured by some observer 
is E = —p-u, where pis the light's 4-momentum and u is the observer's 4- velocity. 
We take this observer to be at rest, so 

E = -pt = -gttp' = p' . (160) 

Put p^ = p^ + Sp^, where p^ describes the light in the absence of a GW, and 
is a GW shift. In the x-arm, p^ = fkv^l, ±1, 0, 0); the signs correspond to before 
and after bounce. We compute the shift 6p'^ using the geodesic equation (fT2|) : 

^ + T^.pp'^p^ = . (161) 

(We use dp'^/dx = in the background to simplify.) Focus on the fi = t compo- 
nent. The only nontrivial connection is T^xx = dth/2. Using in addition the facts 
that p^ = dx^/dx plus j5* = ±p^ reduces the geodesic equation to 

^ = -k%h. (162) 
dt 2^ ^ ^ 

Integrating over a round trip, we find 

6p' = [/l(rround-trip) - HO)] . (163) 

So, we finally find the proper frequency: 

UJproper =E/h = W ^1 + ^ [/l(0) - /l(Ti.ound-trip)]^ 

~ LO . (164) 
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On the second line of Eq. (jl64p . we take Tjound-trip to be much smaher than the 

wave period. This is an excellent approximation for ground-based interferometers; 

the exact result must be used for high-frequency response in space. 

Turn next to the proper round-trip time. The metric (jl53p shows us that 

proper time measured at fixed coordinate is identical to the coordinate time t. 

For light traveling in the x-arm, = —c^dt'^ + (1 + h)dx^ , so 

dx = ±cdt{l-^ +0{h^) . (165) 

Now integrate over x from to L and back, and over t from to Tj-ound-trip: 

^ r Tfound — trip 

'^L = cTround-trip ~ 7^ I 

^ Jo 

c l'"^^/" 

= cr,ound-trip -1, hdt + 0{h^) . (166) 

^ Jo 

We thus find 

_ 2L 1 f^^^\, 

J round— trip — ~^ o / ildt 

C ^ Jo 



_^(i.i„). 

The second line describes a wave which barely changes during a round trip. 
The total phase change is found by combining Eqs. (jl59p . (|164p and (jl67p : 

Ad), = ^ + ^ ^'"^^^ hdt^ (^1 + i [h{0) - h{2L/c)]^ 

2u)L / 1,\ , , 

l + -h] . (168) 



c V 2 , 

The second line is for a slowly varying wave, and the first is exact to order h. We 
will use the slow limit in further calculations. Repeating for the y-arm yields 

= ^ (l - 1/.) . (169) 

Notice that this GW acts antisymmetrically on the arms. By contrast, any laser 
phase noise will be symmetric: because the same laser state is sent into the arms 
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by the beam splitter, we have A$^°^'^'^ = <I>^°'^''. We take advantage of this by 

reading out hght produced by destructive interference at the beamsphtter: 

= 2A$^w . (170) 

An L- shaped interferometer is sensitive only to the GW, not to laser phase noise. 
This is the major reason that this geometry is used; even if an incident wave is 
oriented such that the response of the arms to the GW is not asymmetric, one is 
guaranteed that phase noise wih be cancehed by this configuration. 

From basic principles we turn now to a brief discussion of current and planned 
detectors. Our goal is not an in-depth discussion, so we refer readers interested 



in these details to excellent reviews by 



Hough &: Rowan 



detail the characteristics of the various detectors) and 



20001') fwhich covers in 



Tinto k Dhurandail (|2005l ) 



(which covers the interferometry used for space-based detectors). 
8.2 Existing and planned detectors 

When thinking about GW detectors, a key characteristic is that the frequency 
of peak sensitivity scales inversely with armlength. The ground-based detectors 
currently in operation (and undergoing or about to undergo upgrades) are sensi- 
tive to waves oscillating at 10s - 1000s of Hertz. Planned space-based detectors 
will have sensitivities at much lower frequencies, ranging from 10~^ - 0.1 Hz 
(corresponding to waves with periods of tens of seconds to hours). 

The ground-based detectors currently in operation are LIGO {Laser Interfer- 
ometer Gravitational-wave Observatory), with antennae in Hanford, Washington 
and Livingston, Louisiana; Virgo near Pisa, Italy; and GEO near Hanover, Ger- 
many. The LIGO interferometers each feature 4-kilometer arms, and have a peak 
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sensitivity near 100 Hz. Virgo is similar, with 3-kilometer arms and sensitivity 

comparable to the LIGO detectors. GEO (or GEO600) has 600-meter arms; as 
such, its peak sensitivity is at higher frequencies than LIGO and Virgo. By using 
advanced interferometry techniques, it is able to achieve sensitivity competitive 
with the kilometer-scale instruments. All of these instruments will be upgraded 
over the course of the next few years, installing more powerful lasers, and re- 
ducing the impact of local ground vibrations. The senstivity of LIGO should be 
improved by ro ughly a factor of ten, and the bandwidth increased as well. See 



Fritschell (|2003l ) for detailed discussion. 

There are also plans to build additional kilometer-scale instruments. The de- 
tector AIGO (Australian International Gravitational Observatory) is planned as a 



detec tor very similar to LIGO and Virgo, but in Western Australia (jMcClelland 



2OO2I ). This location, far from the other major GW observatories, has great 



potential to improve the ability of the worldwide GW dete ctor network to deter- 



mine the characteristics of GW events ( Searle et al. 



20061 ). In particular, AIGO 



should be able to break degeneracies in angles that determine a source's sky 
position and polarization, greatly adding to the astronomical value of GW ob- 
servations. The Japanese GW community, building on their experience with the 
300-meter TAMA interferometer, hopes to build a 3-kilometer underground in- 
strument. Dubbed LCGT (Large-scale Cryogenic Gravitational-wave Telescope), 
the underground location takes advantage of the fact that local ground motions 
tend to decay fairly rapidly as we move away from the earth's surface. They plan 
to use cryogenic cooling to reduce noise from thermal vibrations. 

In space, the major project is LISA {Laser Interferometer Space Antenna), 
a 5-million kilometer interferometer under development as a joint NASA-ESA 
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mission. LISA will consist of three spacecraft placed in orbits such their relative 

positions form an equilateral triangle whose centroid lags the earth by 20°, and 

whose plane is inclined to the ecliptic by 60° ; see Fig. [TTl Because the spacecraft 

are free, they do not maintain this constellation precisely; however, the variations 

in armlength occur on a timescale far longer than the periods of their target 



waves, and so can 



} e mod eled out without too much difficulty. The review by 



Tinto &: Dhurandaii (|2005l ) discusses in great detail how one does interferometry 
on such a baseline with time-changing armlengths. LISA is being designed to 
target waves with periods of several hours to several seconds, a particularly rich 
band for signals involving black holes that have 10^ Mq ^ M ^ lO'^M©; the 



LISA Pathfinder, a testbed for soi ne of t 



for launch in the very near future (IVitale 



le in ission's components, is scheduled 



20051 ). 



Somewhat smaller than LISA, The Japanese GW community has proposed 
DECIGO [DECI-hertz Gravitational-wave Observatory), a space antenna to tar- 
get a band at roughly 0.1 Hz. This straddles the peak sensitivities of LISA and 



terrestrial detectors, and may thus act as a. 



one band to the other. See 



Kawamura et al 



jridg e for signals that evolve from 



(|2006l ) for further discussion. 



It's worth nothing that, in addition to the laser interferometers discussed here, 
there have been proposals to measure GWs using at om interferometry. A par - 



Dimopoulos et al. 



(120081). 



ticularly interesting proposal has been developed by 
Sources of noise in such experiments are quite different than in the case of laser 
interferometers, and may usefully complement the existing suite of detectors in 
future applications. 



Compact binary gravitational waves 71 
8.3 Measuring binary signals 

The central principle guiding the measurement of GWs is that their weakness 
requires phase coherent signal measurement. This is similar to how one searches 
for a pulsar in radio or x-ray data. For pulsars, one models the signal as a sinusoid 
with a phenomenological model for frequency evolution: 

$0, /o, /o, /o, . . .) = $0 + 27r (^fot + ^fot^ + ^fot^ + . . }j . (171) 

The cross correlation of a model, cos[<5], with data is maximized when the param- 
eters (<I>Oj /o) /o) /o) • • •) accurately describe a signal's phase. For a signal that is 
N cycles long, it is not too difficult to show that the cross-correlation is enhanced 
by roughly \/iV when the "template" matches the data. 

For binary GWs, one similarly cross-correlates models against data, looking 
for the template which maximizes the correlation. [Given the signal weakness 
implied by Eq. ()152p . the cross-correlation enhancement is sure to be crucial 
for measuring these signals in realistic noise.] Imagine, for example, that the 
rule given by Eq. ([90j) accurately described binary orbits over the band of our 
detectors. We would then expect a model based on 



^t;^c,tc,M) = ^c 



c^{t.-ty"' 



(172) 



bGM 

to give a large correlation when the coalescence phase <l>c, coalescence time to and 
chirp mass M are chosen well. As we have discussed at length, Eq. (I90p is not 
a good model for strong-field binaries. The need to faithfully track what nature 
throws at us has been a major motivation for the developments in perturbation 
theory, pN theory, and numerical relativity discussed here. 

When one determines that some set of parameters maximizes the correlation, 
that set is an estimator for the parameters of the binary. More formally, the cross- 
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correlation defines a likelihood function, which give s the probabiUty of measuring 



some set of parameters from the data (jFinn 



1992I ). By examining how sharply 



the likelihood falls from this maximum, one can estimate how accurately data 
determines parameters. For large cross-correlation (large signal-to-noise ratio), 
this is simply related to how the wave models vary with parameters. Let 9"" 
represent the a-th parameter describing a waveform h. If {h\h) denotes the cross- 
correlation of h with itself, then define the covariance matrix 



dh 
d9^ 



dh_y 



(173) 



(where the —1 power denotes matrix inverse). Diagonal elements of this matrix 
are 1 — a para r neter errors; off-diagonal elements describe parameter correla- 



tions. See iFinnI (jl992l ) for derivations and much more detailed discussion. In the 
discussion that follows, we use Eq. (jl73p to drive the discussion of how model 
waveforms are used to understand how well observations will be able to determine 
the properties of GW-generating systems. 



8.4 What we learn by measuring binary GWs 

8.4.1 Overview Given all that we have discussed, what can we learn by 
observing compact binary merger in GWs? To set the context, consider again 
the "Newtonian" waveform: 

2/3 



2GM 


fuGMf 


c^D 




4GM 


firCMf 


c^D 


[ C3 



C^{tc-t) 

5GM 



j (1 + cos^t)cos2$Ar(t) , 

\ 2/3 

I cost sin 2<I>7v(i) ; 

5/8 



(174) 
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A given interferometer measures a combination of these two polarizations, with 

the weights set by the interferometers antenna response functions: 

h^,,, = F+{e,^,i;)h+ + F,,ie,^,^lJ)h^ . (175) 

The angles {9, (j)) give a sources position on the sky; the angle ip (in combination 
with i) descr ibes the orient ation of the binary's orbital plane with respect to a 



detector. See iThornd (|l987l ). Eqs. (103) - (104) for further discussion 

Imagine that Eq. (jl74p accurately described GWs in nature. By matching 
phase with the data, measurement of the GW would determine t he chirp mass A4 . 



Calcu lations using Eq. (|173p to estimate measurement error (jFinn &: Chernofll 



19931 1 show that ^A should be determined with exquisite accuracy, AAi/^A oc 
1/Ncyc, where A'^cyc is the number of GW cycles measured in our band. 

The amplitude of the signal is determined with an accuracy A A/ A ~ 1/SNR. 
This means that, for a given GW antenna, a combination of the angles 9, (p, t, ip, 
and the source distance D are measured with this precision; however, those pa- 
rameters are not invidually determined. A single interferometer cannot break the 
distance-angle correlations apparent in Eqs. (|174p and p75p . Multiple detectors 
(which will each have their own response functions and F-^) are needed to 
measure these source characteristics. This is one reason that multiple detectors 
around the globe are being built. (For LISA, the constellation's motion around 
the sun makes Fj^ and F^ effectively time varying. The modulation imposed by 
this time variation means that the single LISA antenna can break these degenera- 
cies, provided that a source is sufficiently long-lived for the antenna to complete 
a large fraction of an orbit.) 

As we have discussed, Eq. (jl74p does not give a good description of GWs 
from strong-field binaries. Effects which this "Newtonian gravity plus quadrupole 
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waves" treatment misses come into play. Consider the pN phase function, Eq. 

()127p . Not only does the chirp mass M influence the phase; so too does the 

binary's reduced mass and its "spin-orbit" and "spin-spin" parameters /3 and 

a. The pN phasing thus encodes more detail about the binary's characteris- 

tics. Unfortunately , these parameters may be highly correlated. For example. 



Cutler &: Flanaganl (|1994| ) show that when precession is neglected, errors in a 
binary's reduced mass and spin-orbit parameter are typically 90% or more corre- 
lated with each other. This is because the time dependence of their contributions 
to the phase is not very different. 

These correlations can be broken when we make our models more complete. 
As we've discussed, the spin precession equations ()124p and ()125p cause 3 and a 



to osc illate. This modulates the waves' phase; as first demonstrated 



Lang HughesI (|200g), the mod 



by 



Vecchio 



(|2004l ) and then examined in greater depth by 
ulations break parameter degeneracies and improve our ability to characterize 
the sy stem whose waves we measure. Figure (jl2p . taken from 



Lang Sz Hughes 



(|2006l ). shows this effect for LISA measurements of coalescing binary black holes. 
Accounting for precession improves the accuracy with which the reduced mass is 
measured by roughly two orders of magnitude. We similarly find that the mem- 
bers' spins can be determined with excellent accuracy. GW measurements will 
be able to map the mass and spin distributions of coalescing binaries. 

Equation (jl74p is also deficient in that only the leading quadrupole harmonic 
is included. As the discussion in Sec. 15.2.11 demonstrates, that is .j ust one har 



monic among many 



2007 



that contribute to a bina,ry's G Ws. Recent work (jArun et al. 



Trias &: Sintes 



2008 



Porter &: Cornish 



20081 ) has looked at how our ability 



to characterize a source improves when those "higher harmonics" are included. 
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Typically, one finds that these harmonics improve our ability to determine a bi- 
nary's orientation l. This is largely because each harmonic has a slightly different 
functional dependence on t, so each encodes that information somewhat differ- 
ently than the others. The unique functional dependence of each harmonic on l 
in turn helps break degeneracies between that angle and the source distance D. 
8.4.2 "Bothrodesy" : Mapping BLACK HOLE SPACETIMES Extreme mass 
ratio captures may allow a unique GW measurement: We may use them to "map" 
the spacetimes of large black holes and test how well they satisfy the (rather 
stringent) requirements of GR. As discussed in Sec. [U an extreme mass ratio 
inspiral is essentially a sequence of orbits. Thanks to the mass ratio, the small 
body moves through this sequence slowly, spending a lot of time "close to" any 
orbit in the sequence. Also thanks to the mass ratio, each orbit's properties are 
mostly determined by the larger body. In analogy to geodesy, the mapping of 
earth's gravity with satellite orbits, one can imagine bothrodesy^, the mapping of 
a black hole's gravity by studying the orbits of inspiraling "satellites." 

In more detail, consider first Newtonian gravity. The exterior potential of a 
body of radius R can be expanded in a set of multipole moments: 

^N = + Mi^YUe,^). (176) 

Studying orbits allows us to map the potential <I>Ar, and thus to infer the moments 
Mim- By enforcing Poisson's equation in the interior, V^<I>Ar = AnGp, and then 
matching at the surface R, one can relate the moments M;^ to the distribution 
of matter. In this way, orbits allow us to map in detail the distribution of matter 
in a body like the earth. 



^This name was coined by Sterl Phinney, and comes from the word Po9po<;, which refers to 
a sacrificial pit in ancient Greek. This author offers an apology to speakers of modern Greek. 
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Bothrodesy applies the same basic idea to a black hole. The spacetime of any 

stationary, axisymmetric body can be described by a set of "mass moments" Mi, 

similar to the M/^ of Eq. ()176p : and a set of "current moments" Si which describe 

the distribution of mass-energy's flow. What makes this test powerful is that the 

moments of a black hole take a simple form: for a Kerr black hole (jlOip with 

mass M and spin parameter a, 



Ml + iSi = M{iay . 



(177) 



A black hole has a mass moment Mq = M and a current moment 5*1 = aM 
(i.e., the magnitude of its spin is aM, modulo factors of G and c). Once those 
moments are known, all other moments are fixed if the Kerr solu tion describes 



the spacetime. This is a restatement of the "no hair" theorem (j Carter 



Robinson 



1971 



19751 ) that a black hole's properties are set by its mass and spin. 
The fact that an object's spacetime (and hence orbits in that spacetime) is de- 
termined by its multipoles, and that the Kerr moments take such a simple form, 
suggests a simple consistency test: Develop an algorithm for mapping a large 



object's moments by stud ying or 



moments satisfy Eq. (jl77p . 



)its o f that object, and check that the I > 2 



ment can in principle be done, and 



RvanI (119951) f i rst de monstrated that such a measure- 



Brink 



(|2008l ) has re cently clarified what mus t 



be done for such measurements to be done in practice. 



Colhns fc HughesI tooi ) 



took the first steps in formulating this question as a null experiment (w i th the 



Glampedakis fc BabakI tood ) 



Schwarzschild solution as the null hypothesis), 
formulated a similar approach appropriate to Kerr black holes, and Vigeland 
(Vigeland & Hughes, in preparation) has recently extended the Collins & Hughes 
formalism in that direction. 



A robust test of the Kerr solution is thus a very likely outcome of measuring 
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waves from extreme mass ratio captures. If, however, testing metrics is not your 

cup of tea, precision black hole metrology may be: In the process of mapping a 

spacetime, one meas ures with exquisi t e accu racy both the mass and the spin of 



the large black hole. 



Barack &: Cutleii (|2004l ) have found that in most cases these 



events will allow us to determine both the mass and the spin of the large black 
hole with 0.1% errors are better. 

8.4.3 Binary inspiral as a standard "siren." A particularly exciting 
astronomical application of binary inspiral comes from the fact that the GWs 
depend on, and thus directly encode, distance to a source. Binary inspiral thus 
acts as a standard candle (or "standard siren," so named because it is ofte n usefu. 



to reg ard GWs as soundlike), with GR providing the standardization. 



Schutz 



(|l986l ) first demonstrated the po wer of GW obser vatio ns of merging bin a ries to 



pin down the Hubble constant; 



Markovid (11993) and 



Finn fc ChernofB (11993) 



analyzed Schutz's argument in more detail, in addition assessing how well other 
cosmo logical parameters could be determined. More recently, 



Holz Hughes 



(j2005l ) have examined what can be done if a GW merger is accompanied by an 
"electromagnetic" counterpart of some kind. We now describe how inspiral waves 
can serve as a standard siren. 

Imagine that we measure a nearby source, so that cosmological redshift can be 
neglected. The measured waveform generically has a form 

GM{mi) 



h 



c^r 



-A{t)cos[<^{t;mi,Si)] , 



(178) 



where rrii are the binary's masses, Sj are its spins, and M{mi) is a function of 
the masses with dimension mass. For example, for the Newtonian quadrupole 
waveform (|89p . this function is the chirp mass, M(mj) = M = {mim2)^^^ / {mi + 
m2)^l^ . The function Aii) is a slowly varying, dimensionless function which 
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depends most strongly on parameters such as the source inchnation l. 

Now place this source at a cosmological distance. Careful analysis shows that 



the naive Euclidean distance measure r should be t 



( Carroll 



2004 Chap. 8); see, e.g., 



r e prop er motion distance Dm 



Finn Sz Chernoff 



(|l993l ) for a derivation. Also, 



all timescales which characterize the source will be redshifted: If r is a timescale 
characterizing the source's internal dynamics, r ^ (1 + z)t. 

What is the phase <I> for this cosmological binary? Because it evolves solely due 
to gravity, any parameter describing the binary's dynamics enters as a timescale. 
For example, a mass parameter becomes a time: m — > = Gm/c^ . This time 
suffers cosmological redshift; the mass that we infer by measuring it is likewise 
redshifted: mmeas = (1 + •2)"iiocai- Spin variables pick up a squared reshift factor: 
•^meas = (1 + •2)^5'iocai- This tells US is that redshift ends up degenerate with other 
parameters: A binary with masses rrij and spins Sj at redshift z has a phase 
evolution that looks just like a binary with (1 + z)mi, (1 + z)^Sj in the local 
universe. So, if we put our source at redshift z, Eq. (jl78p becomes 

h = ^^^^^(t) cos [Ht; (1 + z)mi, (1 + zfSi)] . (179) 

Recall that proper motion distance is related to luminosity distance hj Dm = 
Dl/{1 + z). Because we don't measure masses but rather (1 + z) times masses, 
it makes sense to adjust the amplitude and put 

h = + cos [^{t- (1 + z)m„ (1 + zfS,)] . (180) 

The key point here is that measurements directly encode the luminosity distance 
to a source, D^; however, they do not tell us anything about a source's redshift 
z. In this sense GW measurements of merging binaries can be distance probes 
that are highly complementary to most other astronomical distance measures. 



Compact binary gravitational waves 



accuracy using ground-b ased instruments (e.g 



to ^ 1 — 5% from space fiLang fc Hughes 



2008 



Porter &: Cornish 



2009 . 
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d be measured to ^ 10 — 20% 



Cutler & Flanae 


an 


1994 


), anc 


Arun et al. 


2007 




Trias & Sintes 



20081 ). 

Suppose that we measure GWs from a merging compact binary, allowing us to 
measure Dl with this accuracy. If it is possible to measure the source's redshift 
[either f rom th e stati stical properties of the dist ribution of events, as empha- 



sized by 



Schutd (|l986l l and 



Chernoff &: Finn 



an "electromagnetic" event (jHolz &: Hughes 



'll993l ). or by direct association with 



20051 )]. then one may be able to ac- 



curately determine both distance and redshift for that event — a potentially 
powerful constraint on the universe's cosmography with completely different sys- 
tematic properties than other standard candles. An example of an event which 
may constitute such a standard siren is a short-hard gamma-ray burst. Evidence 
has accumulated recently consistent with the hypothesis that at l east some short 
hard bursts are associated wit 



h NS-NS or NS-B: 



Nakar. Gal- Yam Sz Fox 



20061 . 



Perlev et al 



i mergers (e.g., I Fox et al 



2005 



20081 ). Near simultaneous measure- 



ment of a GW signal with a short-hard burst is a perfect example of what can 
be done as these detectors reach maturity and inaugurate GW astronomy. 
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Figure 1: The g eome try of a generic Kerr black hole orbit [taken from 



Drasco &: Hughea (j2006l )]. This orbit is about a black hole with spin parame- 
ter a = 0.998M (recall a < M, so this represents a nearly maximally spinning 
black hole). The range of its radial motion is determined hy p = IGM/c^ {G 
and c are set to 1 in the figure) and e = 1/3; 9 ranges from 60° to 120°. The left 
panel shows the torus in coordinate space this torus occupies. The right panel 



illustrates how a generic orbit ergodically fills this torus. 
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Figure 2: Orbital frequencies for generic Kerr black hole orbits. We vary 
the orbits' semilatus rectum p, but fix eccentricity e = 0.5 and inclination 
parameter ^min = 75°. Our results are plotted as a function of semimajor 
axis A = p/ \/l — e^. All three frequencies asymptote to the Keplerian value 



Q = \J GAi/A^ in the weak field, but differ significantly from each other in the 
strong field. 
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Figure 3: Waveform arising from a generic geodesic black hole orbit, neglecting 
orbital evolution due to backreaction. This orbit has p = SGM/c^ and e = 0.5, 
corresponding to motion in the range 16GM/3c2 <r{t) < 16GM/c^; it also has 
^min = 60°. The large black hole has a spin parameter a = 0.9M. Note that the 
wave has structure at several timescales, corresponding to the three frequencies 
Or, Qg, and (cf. Fig. [2D. 




Figure 4: The ev olution o f circu lar orbits (e = 0) about a black hole with a 



0.8M; taken from [Hughes] (|200d ). The inclination angle l is given by i ~ tt/2 
Oram, the equality is exact for ^min = and for a = 0- In the general case, 



this relation misestimates by ^ 3%; see iHughed (120001 ) for detailed discussion. 
The dotted line is this hole's "last stable orbit"; orbits to the left are unstable 
to small perturbations, those to the right are stable. Each arrow shows how 
radiation tends to evolve an orbit; length indicates how strongly it is driven. 
These orbits are driven to smaller radius and to (very slightly) larger inclination. 
The extremely long arrow at t ~ 120°, r = 7GM/c^ lies very close to the last 
stable orbit. As such, a small push from radiation has a large impact. 
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Figure 5: Plus polarization of wave generated 
massive black hole; this figure is adapted from 



3V a small body spi r ahng into a 



Sundararajan et al. 



([2003). The 



amplitude is scaled to the source's distance D and the small body's mass ^; time 
is measure in units c^/GM. For this calculation, the binary's initial parameters 
are p = lOGM/c^, e = 0.5, and Oram ^ 61°; the binary's mass ratio is fixed to 
fi/M = 0.016, and the larger black hole's spin parameter is a = 0.5M. The insets 
show spans of length At ~ lOOOGM/c^ early and late in the inspiral. Note the 
substantial evolution of the wave's frequencies as the orbit shrinks. 
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Figure 6: The first three contributions to the + GW polarization, and their sum. 
In all panels, we plot {c^D/Gfi)h^ versus c^t/GM. The upper left panel gives 
the sum [cf. Eq. (jl26p ] arising from H^, H^'^ , and H}^; the other panels show 
the individual contributions from those H^. Although subdominant, the terms 
other than make a substantial contribution to the total, especially at the end 
of inspiral (here normalized to c^t/GM = 1). 
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Figure 7: Illustration of precession's impact on a binary's waves. The top panels 
show and /ix for a binary that contains nonspinning black holes; the lower 
panels show the waveforms for a binary with rapid rapidly rotating (a = 0.9M) 
holes. The strong amplitude modulation is readily apparent in this figure. Less 
obvious, but also included, is the frequency modulation that enters through the 
spin-dependent orbital phase parameters (3 and a [cf. Eq. (jl27p ]. 
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Figure 8: Slicing of "spacetime" into "space" and "time." The vector n is normal 
to the slice at t. An observer who moves along n experiences an interval of 
proper time dr = a dt, where a is the lapse function. The shift /? describes the 
translation of spatial coordinates from slice-to-slice as seen by that normal 
observer. 
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Figure 9: Left panel: Comparison between the numerical relativity computed 
frequency and phase and the effective one-body frequency and phase. Right 
panel: Gravitational waveform computed by those two methods. These plots are 
for the coalescence of two non-spinning black holes with a mass ratio 1712/ mi = 
4. Figure kindly prov ided to the author by Alessandra Buonanno, taken from 



Buonanno et al. 



(|2007l ). Note that G = c = 1 in the labels to these figures. 
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Figure 10: Schematic of an interferometer that could be used to detect GWs. 
Though real interferometers arc vastly more complicated, this interferometer 
topology contains enough detail to illustrate the principle by which such mea- 
surements are made. 
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Figure 11: Schematic of the LISA constellation in orbit about the sun. Each 
arm of the triangle is 5 x 10^ km; the centroid of the constellation lags the Earth 
by 20°, and its plane is inclined to the ecliptic by 60°. Note that the spacecraft 
orbit freely; there is no formation flying in the LISA configuration. Instead, 
each spacecraft is in a slightly eccentric, slightly inclined orbit; their individual 
motions preserve the near-equilateral triangle pattern with high accuracy for a 
timescale of decades. 
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Figure 12: Accuracy with which reduced mass fi is measured by LISA for binaries 
at 2; = 1 with masses mi = 3 x 10^ Mq, m2 = 10^ M©. The two curves come from 
a Monte-Carlo simulation in which the sky is populated with 10^ binaries whose 
positions, orientations, and spins have been randomly chosen. Horizontal axis is 
the logarithmic error A/// ji; vertical axis is the number of binaries that fall in an 
error bin. The dashed line neglects spin precession effects; note that the distri- 
bution peaks at an error A/i/iJ, ~ 0.03. The solid line includes spin precession; 
note that the peak error is smaller by roughly two orders of mangnitude. 




